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Abstract 

The Lyman alpha emission line (Lya) of neutral hydrogen (Hi) is intrinsically the brightest emission feature 
in the spectrum of astrophysical nebulae, making it a very attractive observational feature with which to 
survey galaxies. Moreover as an ultraviolet (UV) resonance line, Lya possesses several unique characteristics 
that make it useful to study the properties of the interstellar medium (ISM) and ionizing stellar population at 
all cosmic epochs. In this review I present a summary of Lya observations of galaxies in the nearby universe. 
By UV continuum selection, at the magnitudes reachable with current facilities, only « 5% of the local galaxy 
population shows a Lya equivalent width (Wh ya ) that exceeds 20 A. This fraction increases dramatically at 
higher redshifts, but only in the local universe can we study galaxies in detail and assemble unprecedented 
multi-wavelength datasets. I discuss many local Lya observations, showing that when galaxies show net Lya 
emission, they ubiquitously also produce large-scale halos of scattered Lya, that dominate the integrated 
luminosity. Concerning global measurements, we discuss how W\, ya and the Lya escape fraction (/^“) are 
higher (Whya > 20A and y c a > 10%) in galaxies that represent the less massive and younger end of the 
distribution for local objects. This is connected with various properties, such that Lya-emitting galaxies have 
lower metal abundances (median value of 12 + log(0/H) ~ 8.1) and dust reddening. However, the presence 
of galactic outflows/winds is also vital to Doppler shift the Lya line out of resonance with the atomic gas, 
as high Whya is found only among galaxies with winds faster than ~ 50km s -1 . The empirical evidence is 
then assembled into a coherent picture, and the requirement for star formation driven feedback is discussed 
in the context of an evolutionary sequence where the ISM is accelerated and/or subject to hydrodynamical 
instabilities, which reduce the scattering of Lya. Concluding remarks take the form of perspectives upon 
future developments, and the most pressing questions that can be answered by observation. 

Keywords: Galaxies: evolution - Galaxies: starburst - Physical data and processes: radiative trans¬ 
fer — Ultraviolet: galaxies — Galaxies: individual: Haro 11: iZw 18; SBS 0335-052; Haro 2; Mrk701; 
IRAS 08208+2816; ESO338-IG04. 


1 INTRODUCTION 

The Lyman alpha (Lya) emission line of atomic hy¬ 
drogen (Hi) is intrinsically the most luminous spectral 
emission feature in astrophysical nebulae. It is produced 
by the spontaneous decay from the first excited state 
to the ground state (n = 2 —> 1, where n is the princi¬ 
ple quantum number), where the energy difference be¬ 
tween the levels dictates a photon energy of 10.2 eV, 
or wavelength A = 1215.67 A. After the capture of an 
electron in ionized gas, the transition probabilities of 
the following radiative cascade are such that 68% of re¬ 
combination^ involve the production of a Lya photon 
( Diikstral . 12014 ). Thus if all ionizing photons had an en¬ 
ergy slightly above the ionization edge of Hi, ss 50% of 
the total ionizing energy would be reprocessed into the 
Lya line. Because the line may be intrinsically so lumi¬ 
nous, Lya was initially proposed as the spectral beacon 


by which to identify the first gen erations o f primeval 
galaxies almost five decades ago ([Partridge fe PeeblesL 


1 9671) . 


Lya plays a pivotal role in contemporary astro¬ 
physics, where it is used to identify high-redshift (z) 
star-forming galaxies by either narrowband filter obser¬ 
vations or spectroscopic techniques, and frequently to 
confirm the redshift of candidate galaxies selected by 
other methods. Indeed the recognition that Lya obser¬ 
vations have attained makes the detection of Lya a pri¬ 
mary science goal for many new instruments on large 
telescopes: first light h as rec ently been seen by Hyper- 
Suprime Cam (HSC. iTakada . 2010h at the Subaru Tele- 
scope, the Multi Unit Spectroscopic Explorer (MUSE, 
Bacon et al. iimj at ESO’s Very Large Teles cope, and 


the Cosmic Web Imager (CWI, iMartin et all l2010h at 
Palomar. Furthermore, the Hobby Eberly Dark Energy 
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Experiment (HETDEX, Iffill et all l2008h will find on 
the order of 10 6 Lya galaxies, and Lya detection is 
among the science goals for spectroscopic instruments 
on the James Webb Space Telescope (JWST) and all 
plans for Extremely Large Telescopes (ELTs). 


• The Evolution of Galaxies. The fraction of galax¬ 
ies with Wi,y a > 20A (the canonical definition of 
a Lya-emitter, LAE), at absolute UV magnitudes 
bri ghter than -18, is just 5% in the nearby universe 
( Cowie et al. . l2010i l. where LAEs are rare. However 


onant one, and Lya is also absorbed by Hi in the same 
transition. After absorption to the 2 P level (in the ab¬ 
sence of electron collisions), there is is no alternative but 
for the electron to de-excite through Lya. The optical 
depth of Hi, as seen in the core of the absorption line, is 
given by r 0 w 3 x 10” 14 y/l0 4 /T • IVhi, where T is the 
temperat ure of the gas and Nm the column density 
in cm -2 ( Verhamme et all 1 20061) . Thus at the limiting 
temperatures to which hydrogen can remain neutral, Hi 
becomes optically thick to Lya at IVhi ~ 3 x 10 13 cm -2 . 
Assuming a number density of 1 atom per cubic cm, 
a cloud will exceed r = 1 when its diameter exceeds 
10 -5 pc, or just 2 AU. Taking the Milky Way as an ex¬ 
ample, there are very few sightli nes t hrough which Wn 
drops below 10 2 ° cm -2 ( Kalberla et all 2005 1. imply¬ 
ing that Lya would almost always see upwards of 10 6 
optical depths. 

The upshot is that in most galaxies, Lya undergoes 
a radiative transfer process: photons scatter until they 
either escape from the galaxy or are absorbed by a dust 
grain, and dust extinction is also strongest in the far 
UV. This transfer may be thought of as a diffusion-like 
process, where photons take r andom walks in both phys¬ 
ical and frequency space ( Osterbrocll Il962h . The path 
taken by Lya is entirely regulated by the distribution of 
Hi that it encounters and must traverse, which in turn 
determines the likelihood that Lya will encounter dust. 
Fortunately Lya may see a significantly lower optical 
depth if it is shifted in frequency or the Hi is moving; the 
former can occur either after many scattering events as 
it diffuses in frequency through the redistribution pro¬ 
file, or by scattering in Hi that is itself kinematically 
offset from the Hii media where the Lya formed. Ul¬ 
timately the emitted Lya luminosity (also its EW and 
departure from intrinsic Lya/Ha ratio of 8.7) will be 
a function of Hi distribution, gas kinematics, dust con¬ 
tent, and galaxy viewing angle. 


1.1 Key Applications of Lyman alpha 

Lya transfer makes it hard to interpret the observed 
Lya flux and EW from an individual galaxy, because 
the escape fraction, flff . a , is difficult to predict for given 
configuration. However the transfer process, and the 
sensitivity of Lya to different ISM properties, is also 
one of the major advantages of the transition. I now 
outline some key applications. 


tion are incomplete. While Lya forms in astrophysical 

to ~ 25% at * r 

i 3 (ShaDlev et al. 

20031) to over 50% 

nebulae at the stated intensities, the transition is a res- 

at z = 6 — 6.5 

Stark et al., 201C 

: Curtis-Lake et al.. 


2012h . fg/ c a ev olves even m o re strongly over the sam e 


redshift range (IHaves et all 1201 lb iBlanc et all 120111) . 


This monotonic evolution, that spans a factor of 100 
in flff a , is a key result of many Lya surveys but has 
no conclusive explanation. Dust and Hi covering have 
both been suggested, and the answer must indeed lie 
among the quantities mentioned above, or combinations 
thereof, in the co-evolving properties of stars, gas and 
dust. 

• The Epoch of Reionization. Lya emission offers 
a unique opportunity to study the ionization state 
of intergalactic medium (IGM) at an epoch where 
other methods - e.g. the Gunn-Peterson trough in 
QSO spectra and Thompson scattering of the cos¬ 
mic microwave background (CMB) - are insensitive. 
The Lya emitter fraction and fljf. a evolution discussed 
above reverses a f ter z ~ 6.5, decreasing to < 20% (e.g. 
Pentericci et all 2014|) at z~ 7. Possible interpreta¬ 


tions inclu de an i ncrease i n the ion izing photon escape 
fraction ( Diikstra fe Jeeson-Daniel . 20131 which is any¬ 
way needed for reionization), but also that and increas¬ 
ingly neutral IGM starts to absorb the Lya produced 
by the galaxies themselves. Disentangling the scenarios 
requires more information and solid constraints on the 
reionization history and topolog y require much larger 
samples (e.g. Jensen et ah, 2014), but these will become 
available in the coming years with HSC. Moreover, as 
the IGM becomes neutral, the damping wing of Hi Lya 
absorption may begin to affect the profile shape of the 
Lya line that is transmitted, giving Lya another unique 
application. 

• Galaxy Kinematics. Lya scatters coherently in the 
restframe of the Hi atom, and at scattering events is 
shifted in frequency by the velocity of the scattering 
medium. Thus as Lya may escape from galaxies be¬ 
cause of frequency shifts, the kinematic structure of the 
atomic gas becomes imprinted onto the line. This man¬ 
ifests as both a redshift (for outflowi ng gas) of the cen¬ 


troid of the main emission peak (e.g. lHashimoto et al 


2013h and also as characteristi c features that modify th e 


shape of the line profile (e.g. Verhamme et al. . 2008h . 


While there are many probes of kinematics in astro¬ 
physics, nebular line kinematics exclusively traces the 
warm ionized medium. Lya kinematics on the other 
hand is shaped by kinematic differences between Hii 
(production) and Hi (scattering) media. Moreover Lya 
is intrinsically very bright, and can be seen redshifted 
from the most distant galaxies. This again provides 
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unique insights into the evolution of the ISM of galaxies, 
which can only feasibly be done with Lya. 

• Atomic Gas Surrounding Galaxies. As well as mod¬ 
ifying the line profile, scattering also chan ges the sur¬ 


face brightness pr o file o f emitted Lya (|Steidel et al 


l201lk lHaves et al. . 2013h . The mechanisms by which 
galaxies obtain the gas they need to fuel star forma¬ 
tion is one of t he most pre s sing issues in extragalac - 
tic science (e.g. Keres et all . 12005 : Dekel et al. . 2009), 
and necessitates a knowledge of the Hi distribution 
outside of star-forming regions and into the circurn- 
galactic medium (CGM). Probing this circumgalac- 
tic Hi is observationally very challenging. In principle 
it can be done by A = 21 cm observations of Hi di¬ 
rectly, but current telescopes cannot push such tech¬ 
niques beyond the very local universe. An alternative 
is to use absorption spectroscopy of background QSOs 
that pierce galaxy halos at different impact factors, en¬ 
abling us to m easure Hi temperatures, densities and 
kinematics (e.g. Lanzetta et allIl995t iTumlinson et al 


120131 : iDanforth et al. . 2014h . Unfortunately appropri¬ 
ately bright QSOs are rare and thus studies, while rich 
with information, are limited to statistical studies of 
the average galaxy. A promising third method is to il¬ 
luminate the circumgalactic Hi with Lya produced in 
the central star-forming regions. Indeed Lya is perfect 
for such an application, being both the brightest intrin¬ 
sic emission line, and being resonant in precisely the 
medium we need in order to image the CGM. 


1.2 Lyman alpha Observations of the Nearby 
Universe This Review 

Section [IT] presents the main astrophysical applications 
of Lya emission, both as a diagnostic of the galaxies 
themselves and the IGM. The key difficulties of ob¬ 
serving high -2 galaxies are that fluxes are low and high 
signal-to-noise data are hard to obtain, that spatial in¬ 
formation is minimal or absent, and that important fea¬ 
tures are redshifted away from atmospheric transmis¬ 
sion windows. In the local universe surface brightness is 
higher by a factor of (1 + 2 ) , and spatial sampling can 
become almost arbitrarily high. Moreover, only in the 
local universe can we assemble the complete set of multi¬ 
wavelength observations, including but not limited to, 
all the continuum bands that probe both hot and cold 
stars, emission lines that provide a wealth of intrinsic 
diagnostics of the nebulae in which Lya forms, direct 
measurements of far infrared continuum for both hot 
and cold dust, direct Hi measurements at A = 21 cm, 
X-ray observations of coronal gas, and many more. 

Indeed the science objectives discussed in Section [TTI] 
can, with the exception of reionization, all be under¬ 
taken in the low -2 universe. Here Lya provides a unique 
suite of information about the ISM of galaxies that 
still cannot be extracted using other techniques. This 
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makes Lya an import observable to obtain in any thor¬ 
ough study of (particularly star-forming) local galaxies. 
Moreover the question can be inverted: when as much 
information on the dust and gas content (distribution, 
kinematics, etc) has been assembled, observations of 
Lya can then be used to calibrate our understanding of 
the Lya transport mechanisms, and the effects of dust, 
gas, star-formation evolutionary stage may all be disen¬ 
tangled. This is the way, for example, we will assemble 
the relevant knowledge to interpret the evolution of the 
Lya fraction with redshift. In turn we will be able to 
calibrate Lya for high -2 galaxy surveys, which will soon 
deliver ~ 1 million objects, by using local Lya emitters 
as analogues - laboratories in which to dissect in de¬ 
tail the processes ongoing in high -2 systems. Local Lya 
observations will allow us to read kinematic informa¬ 
tion off the line profile and conversely to predict the 
flux, EW, and line profile shapes, precisely as needed 
to address topics such as cosmic reionization that really 
hinge upon knowing the spectral profile. Indeed this is 
one of the major legacies established by our ultraviolet 
satellites; the only difficulties are that such satellites are 
both expensive and competitive. 

This review focuses mainly upon empirical studies of 
Lya emission and absorption in star-forming galaxies in 
the local universe. Somewhat arbitrarily, I have defined 
the ‘local universe’ to mean redshifts where space-based 
facilities are needed to observe Lya. In principle this 
means 2 < 1.7 or so, but the most distant samples dis¬ 
cussed are at 2 ~ 1, and thus we are considering roughly 
the latter half of cosmic time. Where appropriate I may 
concentrate upon what about galaxies teach us about 
Lya, or about what Lya teaches us about galaxies. The 
layout of the remainder is as follows: 

• In Section [3] I present a brief history of Lya observa¬ 
tions in the local universe, which were ongoing at a time 
when the first generations of high -2 searches were also 
beginning. This concerns the first vacuum UV observa¬ 
tions of active galactic nuclei (AGN) and star-forming 
galaxies using low dispersion spectrographs. 

• In Section []] we discuss how the Hubble Space 
Telescope changed the landscape by providing high- 
resolution spectra that can resolve the Lya feature and 
interstellar absorption lines, thereby probing atomic gas 
kinematics and covering. 

• Section [5] is concerned with Lya imaging observations, 
also from HST, that simultaneously resolve very fine 
structures and reveal large-scale, diffuse Lya halos. 

• In Section EDI present a large number of key results 
from survey data, that aim to answer questions about 
how various globally measured properties influence Lya 
emission and under what conditions Lya can be ex¬ 
pected to be bright. 

These Sections ED ED aim to establish empirically how we 
have arrived at the current state-of-the-art. 

• In Section [7] I then synthesize all the observational 
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data from the previous Sections, and introduce some 
more speculative discussion about how various pro¬ 
cesses fit together. 

• Finally I do not present explicit conclusions, but close 
the review in Section [8] with a number of perspectives 
and pressing open questions. These are concerned fu¬ 
ture observations and uses of Lya at both low- and 
intermediate-z, with a view to understanding galaxy 
formation. 


2 LYMAN ALPHA OBSERVABLES 

The literature makes use of several observable quanti¬ 
ties that pertain to Lya, that often are specific to how 
much escapes or aim to illustrate how ‘strong’ the fea¬ 
ture is. This Section summarizes some of these quanti¬ 
ties and conventions, and lists several caveats that may 
be considered while reading. 


2.1 Flux, Luminosity, and Equivalent Width 


Lya is mainly produced by recombinations in photoion- 
ized nebulae, where under standard Case B assumptions 
68% of ionizing photons absorbed by hydrogen are re- 
proc essed int o Lya in the following radiative cascade 
(See iDiikstraJ 12014 for an intuitive explanation). For 
continuously star-forming galaxies, with constant star 
formation rate (SFR), the Lya equivalent width (EW) 
ranges between « 80 A. For very young syst ems the EW 
peaks around 300A ( Chariot fe Fall . 1993 1. while very 
high EWs that exceed 1000 A may in principle be ex¬ 
pected f or very low metallicity, population ill stellar 
systems ( SchaereA 1 200. 'it iRaiter et ah . 2010h . Naturally 
if the SFR is declining the intrinsic Lya EW may take 
any value less than these, and it is worth noting that for 
a simple stellar population (SSP), W-[, vn exceeds 20 A 
only during the first 6 Myr ( Leitherer et al. . 1999h . 

Measuring the flux and EW of most emission lines 
is straightforward. However for Lya this is not neces¬ 
sarily so and these quantities may depend upon both 
methodology and definition. As a resonant transition, 
both nebular Lya and continuum radiation with wave¬ 
length A = 1216 A may be absorbed by Hi. Depend¬ 
ing upon the column density, this absorption may reach 
equivalent widths of several tens of A, which is a sub¬ 
stantial fraction of the nebular flux. Fluxes measured 
in a given aperture may or may not be reduced by this 
amount. Narrowband imaging, for example, needs to be 
continuum-subtracted and therefore measures the sum 
of nebular emission and absorption. Spectroscopic ob¬ 
servations, on the other hand, may enable the observer 
to isolate the components and separate nebular emis¬ 
sion from ISM absorption, should this be the goal of 
the measurement. However even with in spectroscopic 
mode isolating the emission-only flux depends upon the 
spectral resolution, and the separation will be much 


easier with high-resolution slit spectrographs than low- 
resolution survey telescopes. 

Not only is Lya absorbed in the ISM but also, de¬ 
pending on temperature and the properties of their 
winds, in the atmospheres of stars. For the hottest O 
stars Lya EWs may be negligible, but as the population 
ages or the star formation rate declines, the dominant 
sourc e of UV continuum wi ll shift to progressively cooler 
stars. Valls-Gabaudl ( 1993ll showed that Lya measure¬ 
ments from some local galaxies may be subject to signif¬ 
icant uncer tainties from ste ll ar absorpt ion, and recent 
models by [Pena-Guerrero fe Leitherer] ( 2013t ) demon¬ 
strate that stellar Lya absorption may reach EWs of 
—30 A. For example, the effect of stellar absorption 
may also vary from O-star dominated systems where 
the nebular EW is high and the stellar feature is negli¬ 
gible, to later B-star systems where the reverse is true. 
The stellar feature may therefore shorten the timescale 
over which an episode of star-formation remains bright 
in Lya, and the effect may even be seen on resolved 
scales within a galaxy. 


2.2 Lya/Balmer Ratios and Escape Fraction 

Equivalent widths have the advantage that only a short 
range in wavelength is needed to make the measure¬ 
ment, over which the effects of interstellar dust (redden¬ 
ing curve, total extinction) should have a negligible ef¬ 
fect. As discussed above, the evolutionary phase of star 
formation dominates the intrinsic EW. Lya/Balmer line 
ratios, however, may also be used as a measure of the 
strength of Lya, and because the intrinsic ratios are 
limited to a narrow range of values, have other ad¬ 
vantages. For example the Ha line (A = 6564. 6lA) is a 
well-known, calibrated tracer of current sta r-formation 
activity in nearby galaxies (e.g. iKennicutt . Il983l ). For 
Case B nebulae at temperatures in the range 5,000- 
20,000 K and electron density in the range n e = 10 2 - 
= 10 4 cm~ 3 , the Lya/Ha ratio ranges between 8.1 and 
11.6 ( Hummer fe Storevl . Il987l) . Thus deviations from 
the intrinsic line ratios encode information about the Hi 
scattering and dust absorption. While different authors 
do tend to adopt slightly different values, the range 
of permitted values is relatively narrow. For this re¬ 
view we will adopt the value of Lya/Ha= 8.7, which 
for T = 10 4 K gas corresponds to n e 


350 cm A and 


as a c onvention can be traced back to at least IHu et al 
(Il998h . 


Frequently we make reference to the Lya escape frac¬ 
tion (/ e L s r)> in an effort to find a quantity that most 
closely reflects the combined impact that gas and dust 
have on suppressing the emission line. We define /^“ as 
the ratio of the emitted Lya luminosity to that intrinsi¬ 
cally produced, but naturally a number of assumptions 
can enter our estimates of the intrinsic value. The most 
robust estimates of fg/ c a will naturally come from com- 
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paring Lya with other hydrogen recombination lines, 
where in practice Ha is the strongest and easiest to 
observe. Assuming that Ha is unobscured, /^ a will 
simply be the observed Lya/Ha ratio divided by its in¬ 
trinsic case B value (8.7 as mentioned above). Of course 
Ha can also be significantly absorbed by dust, and in 
local ‘normal’ galaxie s suffers about 1 magnitud e of ex¬ 
tinction on average ( Kennicutt fe Kentl . Il983 ). Obvi¬ 
ously redder hydrogen lines would be better as they 
suffer less extinction but also become systematically 
weaker in flux, and become harder to observe in the 
near infrared. Radio recombination lines would be op¬ 
timal, suffering no extinction at all, but are even more 
challenging to observe beyond the very local universe. 
Thus often the best route to is to dust-correct Ha 
using the H/3 line, which should recover all the star for- 


mation down to moderate optical depths (IHaves et al 


20051 : lAtek et all |2009aj). However when Ha becomes 


very optical ly thick, in very dusty star-forming galaxies 
(e.g.[Martin ct al.. 120151) . even dust-corrected Ha traces 
only a small fraction of the true ionized gas, making the 
inferred escape fraction prone to strong biases. In such 
instances, one may do better by comparing the Lya- 
derived SFR with that estimated from dust emission in 
the FIR, under the assumption that systematic errors 
on the SFR calibrations are smaller than the fraction 
of Ha that is recoverable. In the highest optical depth 
regimes this is probably true. 


3 THE FIRST LOW REDSHIFT LYMAN 
ALPHA OBSERVATIONS 

The major contributing observatories in this field are 
just three: the International Ultraviolet Explorer (IUE), 
the Hubble Space Telescope (HST), and the Galaxy Evo¬ 
lution Explorer (GALEX), while piecemeal observations 
have been contributed by other telescopes. Restframe 
FUV observations began with the IUE in 1977, HST 
started operation just 13 years later and is still return¬ 
ing Lya data at z < 1 some 25 years on; the demarca¬ 
tion between history and present is arbitrary of course. 
For the sake of this review, I will adopt the pre-HST 
era, which is almost exclusively the IUE and thus is 
also restricted by method to spectrophotometry of pre¬ 
selected targets. In the ‘modern’ era we then have two 
main operational tools. Firstly HST, which like IUE 
performs pointed observations of individual targets, but 
in doing so yields data that are always rich with features 
since both spatial and spectral resolutions are signifi¬ 
cantly higher. Secondly, the GALEX satellite is similar 
to the IUE in resolving power but its strength comes 
instead from its 1.2 degree field of view, which provided 
the survey efficiency to yield the statistical significance 
that was not available at low -2 with any other telescope. 


3.1 The First Lya: Spectra of Active Galactic 
Nuclei 

Even before IUE some extragalactic objects were ob¬ 
served with sounding rocket experiments, which were 
able to launch small UV telescopes above enough of the 
atmosphere to observe in the UV. Specifically target¬ 
ing the first known quasar 3C 273, rocket payloads pro¬ 
vided the first measurement of Lya and B aimer emis- 


sion lines from any astrophysical body (Davidsen et al 


1977 ). Combining the measured Lya flux with optical 
measurements revealed a Lya/H/3 ratio of 4, a mea¬ 
surement that then sat in stark contrast to the value of 
ss 40 that was expected for nebular recombinations and 
a Lya enhancement from collisions. In 1977 Davidsen 
et al. discussed the order-of-magnitude Lya deficit with 
the following statement: 


The most obvious explanation, that the ultraviolet 
lines are attenuated by absorption by dust similar to 
that observed in the interstellar medium, seems un¬ 
tenable in view of observations of Paschen a that in¬ 
dicate the hydrogen lines are unreddened. But, dust 
which is distributed within the line-emitting gas might 
destroy La without having much affect on the Balmer 
and Paschen lines if the nebula has high optical depth 
to La photons. 

[...] 

Whatever mechanism is at work, an understanding 
of the reduced Lya/Hf3 ratio may lead to a vastly im¬ 
proved knowledge of the physical conditions in QSO 
envelopes. 


While this review is not focused upon QSOs, replacing 
“QSO” with “galaxy” in this statement comes precisely 
to the point. Today Lya observations play a major role 
in understanding the interstellar and circumgalactic me¬ 
dia of galaxies at effectively all redshifts. 

The 1978 launch of IUE opened up the restframe 
UV to systematic study. IUE had a single 20"xl0" en¬ 
trance aperture and Short and Long- Wavelength Prime 
channels (SWP and LWP, respectively) that could pro¬ 
vide R ~ 250 spectroscopy between 1150 and 3000 A. 
First observations were again turned to QSOs and var¬ 
ious Seyfert galaxies, and immediately showed Lya 
to be weaker than expected for recombination theory, 
with Lya/H/I and Ha/H/3 values almost never falling 
along characteristic reddening cur ves, and L ya a l ways 
being preferent i ally suppressed (Okefc Zimmermanl . 
Il979l : IWu et al' . 19801 : lLacv et, ah . 1982h . Interoreta- 
tions varied, with suggestions that multiple scatterings 
of Lya could be responsible, that broad-line regions ex¬ 
perience a wide variation in their extinction laws, ion¬ 
ization/excitation from already-excited states, and that 
there may be no re presentative intr insic spe ctral shape 
for Seyfert galaxies (jWu et allll980lh Indeed Lacv et al 
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( 19821 ) concluded that the low observed Lya fluxes were 
likely the result of a combination of reddening, high 
densities, and high Hi optical depth. All these effects 
conspire in the same direction but as no single domi¬ 
nant quantities could be identified, it started to become 
clear that Lya transfer in true astrophysical objects is 
a complicated multi-parametric process. 


3.2 Star-forming Galaxies 


When IUE was first turned to star-forming galaxies, the 
results directly mirrored those obtained in both nearby 
QSOs and also those being reported from hig h -z blind 
narro wband and spectroscopic surveys (see iPritcheti . 
1994 for a review, and the forthcoming PASA review 
in this series by S. Malhotra): Lya was either absent or 
unexpectedly weak in all galaxies. 

With the intent of studying the ana l ogues of primeval 
galaxies at high- 2 , Meier fe TerleviclJ ( 198lh found Lya 
in emission from just one of three nearby Hil-selected 
galaxies, and in that single case with a flux well below 
that expected for the nebular dust content (Figured]). 
This led them to conclude that under normal condi¬ 
tions, Lya emission would be an unlikely phenomenon. 
It was determined that if a normal prescription for dust 
attenuation were used to explain the Lya/H/3 ratio, 
this would greatly over-predict the Ha/H/3 ratio com- 
pa red to observation . Simi lar conclusions were reached 
by lHartmann et al.l ( 19841 ), who furthered the discus¬ 
sion, showing how a mixed medium of Hi and dust could 
preferentially suppress Lya, and that static Hi columns 
of density above 10 19 cm -2 would be needed to recon¬ 
cile the line ratios. Further, they raised the issue that, 
where Hi 21 cm data are available, most blue compact 
galaxies (BCGs) are found to sit inside extended Hi ha¬ 
los of sufficient column density to reproduce the mea¬ 
sured fluxes. All signs pointed towards the fact Lya 
emission would not be the good observational beacon 
to identify pri meva l galaxies in the early universe that 
Partridge &J/eebles (196^) had predicted. 

While the influence of Hi on Lya visibility was 
starting to be seen empirically, some correlation with 
the dust abundance should still be expected, albeit 
with a large spread. After subsequent data acquisi¬ 
tion, the anticorrelation between L ya EW and gas- 
phase metallicity (Z) was discovered (lHartmann et aT 


1988: Calzetti fc Kinnevl Il992t : iTerlevich et al. . 1993; 


Chari ot &l Fall!. 119931) , seemingly confirming the predic¬ 
tion. However the spread on the relation was worryingly 
large, and it was clear tha t something beyond v ariations 
in the extinction law (e.g. Valls-Gabaud . 1993h were be¬ 
hind the wide range of line ratios, with ISM geometry 
a nd holes being th e most often cited scenarios. 

Giavalisco et al. ( 1996] ) performed a full reanalysis 


of the IUE archival data, resulting in several changes. 
Firstly the IUE data reduction software reached a 
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Figure 1.: The first Lya spectra of Hu galaxies. Main 
transitions and cosmic ray hits are labelled. The upper 
panel shows Mrk701 (C0842+163), with a redshift of 
z = 0.0522 and a metallicity of Z « 0AZ@. At the wave¬ 
length of A = 1278A, the galaxy shows only a Lya ab¬ 
sorption feature. The lower panel shows a very different 
object - C 1543+0907, with a redshift of z = 0.0366 and 
metallicity of 1/10 the solar value. At the wavelength of 
A = 1260A, the bright high-contrast Lya emission line 
is obvious. Reproduced by permission of the AAS, from 
Meier & TerlevicH \l98i) . 


higher level of maturity, and new spectral extraction 
and cosmic-ray removal tools were implemented: this re¬ 
duced the measured W^ ya significantly in about half the 
sample, and completely removed the weak Lya feature 
that had been reported in some galaxies. Secondly, these 
authors performed proper spatial matching between the 
IUE and apertures used for optical line spectroscopy. 
When the homogenized reanalysis was complete, both 
the Wi, ya and Lya/Balmer line ratios showed no signif¬ 
icant correlation with nebular dust attenuation or the 
UV continuum colour (/3), and only a weak but signif¬ 
icant correlation with nebular oxygen abundance. Re¬ 
sults concerning the Lya/Ha line ratios can be seen 
Figure [2] 

Regarding these correlations, it is not clear why Lya 
throughput should be more strongly correlated with 
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Figure 2.: The observed ratio of Lya/Ha fluxes, shown on a linear scale against nebular metallicity an d the Ha/H/3 
line ratio for the full sample of star-forming galaxies observed with the IUE, and reprocessed bv lGiavalisco et aT 


( 199fil h Negative Lyq/H a ratios are the res ult of the Lya feature being d o minat ed by I SM absorption. The ref ¬ 
erenc e coding is MT81=lMeier fc r hrlevichl (1 198 lib H 84= Hartmann et alJ ( 1984 ), H88= Hartman n et alJ (1 198 
C94= Calzetti et all ( 1994lh CK92= Calzetti &; Kinnevl (Il992h . D86=lDeharveng et ahl ( 1986h . T93= Terlevich et al 


(1993). The dotted black line shows the effect of the Calzetti et al. ( 2000h dust attenuation curve, assuming intrinsic 
ratios of Lya/Ha=8.7 and Ha/H/1=2.86 and no scattering. Most of the galaxies lie far below this curve. 


oxygen abundance than nebular attenuation. Metals 
alone do not absorb Lya, which can only happen by 
interaction with dust grains. Thus if the W\jy a ~Z anti¬ 
correlation results from a positive correlation between 
metal and dust abundance, then a tighter correlation 
between Lya/Ha and Eb~v would be expected. This 
correlation is completely absent. Furthermore, when the 
dust reddening measured from the B aimer decrement is 
used to correct the observed Lya for extinction, the Lya 
flux does not reach the expected case B recombination 
value in a single galaxy in the IUE sample. This demon¬ 
strates that either some preferential attenuation of Lya 
must be at play in every galaxy, or/and that dust at¬ 
tenuation laws, when applied as a screen of absorbing 
material, are not representative. 

When interpreting these analyses of the IUE samples, 
it is important to keep in mind the selection functions 
by which the individual samples were established. When 
effective telescope areas are small and exposure times 
need to be long, the result is small samples. Some of 
the blue compact dwarf (BCD) and Hii galaxy stud¬ 
ies were designed to study the analogues to ‘primeval’ 
galaxies that are undergoing their first phase of star- 
formation, and the samples d id not include galaxie s 
with strong starbursting nuclei (Hartmann et al.. ll984lh 
Yet the primeval stages of galaxy formation, prior to the 
initial dust and metal production, are expected to be 
short-lived because the first generations of supernovae 


will enrich the local ISM on timescales of just a few Myr. 
Consequently the bulk of the galaxy population we can 
observe in the high redshift universe is likely to be sig¬ 
nificantly m ore metal-enriched than that of p r imeva l 
galaxies (e.g. IPettini et, all 120021 : Ishanlev et all 2003 b 
unless observations catch galaxies in very narrow time 
window. Thus while providing very interesting astro- 
physical laboratories, the samples are biased and not 
necessarily in a direction tuned to the real importance 
of Lya: probing the faint population of normal galaxies 
at z > 2. The IUE samples contain few galaxies that 
can be considered the local analogues of high- 2 ; Ly¬ 
man Break Galaxies (LBG), Lya-emitters, or primeval 
galaxy building blocks. 


4 HIGH RESOLUTION SPECTROSCOPY 
WITH THE HUBBLE SPACE 
TELESCOPE 

High-resolution ultraviolet observations in the nearby 
universe, both imaging and spectroscopy, are one of the 
major legacies of HST. The main UV spectrographs on 
HST were designed to have high resolving powers, which 
necessitated either small entrance apertures or the use 
of narrow slits. The natural consequence is that while 
HST spectra may be rich with features, any measure¬ 
ments strictly reflect the properties of the stars and gas 
that fall within the aperture. In most cases these will 
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Figure 3.: The Lya spectra of first d warf g ala xies t o be o bserved in the far UV with HST. The left panel shows 
iZwl8, first observed with GHRS by Kunth et aJj (119941) . The geocoronal Lya line has been masked out and is 
shown by a line set to zero flux at velocities of —500 to —2000 kms -1 . No intrinsic L ya emissi on i s seen , and 
damping wings are visible that extend to at least 6000 kms -1 . Figure is taken from iMas-Hesse et al. (2003). The 
right panel shows Haro 2, first observed by Lecnieux et al, 019951 ). Haro 2 has a much higher dust abundance than 
I Zw 18, but shows a strong Lya emission line. Furthermore the line profile is P Cygni-likc, which indicates the 
intrinsically p roduced frequencies have been redistributed by scattering in an expanding neutral medium. Data are 
replotted from iMas-Hesse et al.1 (2003). 


be the regions of highest UV surface brightness, which 
in general will be unobscured massive stellar clusters. 
However it is not necessary that these local properties, 
whether regarding Lya or interstellar gas, are represen¬ 
tative of the entire galaxy. Thus we must always keep 
in mind that HST derived measurements of gas kine¬ 
matics, densities, and covering, such as those discussed 
in this Section, are local. 

The first-generation instrument, the Goddard High 
Resolution Spectrograph (GHRS) increased the resolv¬ 
ing power over that of the IUE (R ~ 250) by a factor of 
10-100. This enabled studies of the kinematics and cov¬ 
ering of the atomic gas which are topics completely ab¬ 
sent from discussion in every previous paper presenting 
IUE data. Furthermore, with the IUE observers needed 
to target galaxies at high enough radial velocity to sep¬ 
arate the intrinsic Lya feature from the bright Lya line 
produced in the Earth’s corona, which is brighter than 
any known extragalactic object. Indeed because of the 


tentative anticorre l ation between WLya and metallicity, 
Meier fe Terlevichl ( 198ll) already commented upon how 


it was unfortunate that the lowest metallicity galaxy 
known - lZwl8, with 12 + log(0/H) ss 7.2 - had at too 
small a redshift to be observed with IUE, for logically it 
must be very bright in Lya. With HST/GHRS it could 
finally be observed. 


4.1 HST Finds Deep Lya; Absorption 


GHRS observations of I Zw 18 revealed a profile showing 
only damped Lya a bsorption and no hint of emission 
( Kunth et al. . 19941) . The Lorentzian wings of the ab¬ 
sorption p rofile can be traced out to at least 6000 km s -1 
(Figure El IMas-Hesse et alJfeooJ ), implying a Hi column 
density above several 10 21 cm -2 . Furthermore, mea¬ 


surements of the low ionization stage (LIS) metal lines 
that form in the neutral ISM (e.g. Oi A1302, Sin A1304) 
show that the Hi gas is static with respect to to the Hil 
regions, and at the measured column density Lya ra¬ 
diation at line-centre will have to traverse 10' optical 
depths in order to escape if the Hi is homogeneous. 

Subsequent observations of I Zw 18 with the Space 
Telescope Imaging Spectrograph (STIS) enabled 
detailed, spatially resolved, and empirically well- 
constrained studies of Lya radiative transport, which 
indeed shows that this profile can be reproduced, 
including spatial variation in the damping wings, 
using column densities of Aht ~ 3 x 10 21 cm -2 , and 
Eb-v— 0-05 ( Atek et all 2009bi ). This is fully consis¬ 
tent with the directly observed values on both Nui 
and Eb-v- Moreover the same transport simulations 
predict that for static gas with this Hi column density, 
almost all the Lya radiation is absorbed by the small 
amount of available dust (escape fraction the order of 
10 -4 to 10 -3 ), as many scattering events increase the 
probability of absorption. 

With the HST servicing mission 4, the Cosmic Ori¬ 
gins Spectrograph (COS) was installed on HST. COS 
has a similarly sized entrance window to GHRS and 
similar grating specifications, but can observe much 
larger wavelength range in a single observation, and is 
many times more sensitive. Recently obtained COS ob¬ 
servations of I Zw 18 enable us to go much deeper than 
was previously possible, and a COS/G130M exposure 
of 29 ksec actually does reveal a small bump of Lya 
emission, hidden a t the bottom of the absorption profile 
(Figu re IH see also iLebouteiller et al. 2013 ; James et al.1 


2014I) . The Lya emission feature lies in the centre of 


the absorption trough, and the flux amounts to just 
w4x 10 -16 erg s -1 cm -2 . Recalling at this point that 
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Figure 4.: Very deep Lya spectrum iZwl8 obtained 
with the Cosmic Origins Spectrograph. While previ¬ 
ously referred to as a Lya absorbing galaxy, I Zw 18 
shows a very small but significant bump of Lya emis¬ 
sion. This bump is redshifted from the expected sys¬ 
temic velocity (the dot-dashed vertical line) derived 
from the Ha line. The overlaid red line shows the profile 
of Lya transferring through a static shell of atomic hy¬ 
drogen with a columndensity_of log[lVHi/cm _2 ) = 21.1, 
taken from the ISchaerer et al. ( 201 lh grid of transfer 
models. The strong emission feature is the geocoronal 
Lya emission line. The inset shows the wider spectral 
profile. 


HST spectrographs measure only very local properties, 
we would need to make a large aperture correction to 
estimate the total Lya output. However assuming an ex¬ 
ponential Lya surface brightness profile we would need 
a scale length of 80 arcsec in order for 100% of the Lya 
to be emitted (integrating to infinity). While this is not 
ruled out - Hi is ext ends over at least a square arcmin 
( van Zee et all Il998l ) - this scale length is 40 times the 
UV effective radius, which is an extreme extension of 
Lya compared with other nearby objects (Section [5]). 
The alternative is that a large fraction of the Lya pho¬ 
tons are absorbed by dust after numerou s scattering 
events, as suggested by Atek et al. ( 2009b ). 

In I Zw 18 the Lya bump is also offset from the sys¬ 
temic velocity (measured from Ha) by 350 kms -1 while 
the neutral gas shows bulk velocities of —10 kms -1 
(measured from Sin). The bulk motion is insufficient 
by an order of magnitude to kinematically redshift the 
Lya by such a velocity, so other physical processes must 
be at play. The red line in Figure 0] sho ws a radiat ive 
transp ort model produced by the grid of Schaerer ct al 


( 20111) . for a completely static and dust-free Hi shell 
with a column density log(iVHi/cm _2 ) = 21.1, which 
is within 0.25 dex of previous estimates based upon 
modeling just the absorption feature. A plausible in¬ 


terpretation is that the bump is the red half of a 
double-peaked profile, which results from wing scatter¬ 
ing events that shift photons in frequency many Doppler 
widths into the redistribution profile. If so, and the scat¬ 
tering medium is completely static, there would be a 
corresponding blue peak at -350 kms , this however 
would be hidden below the geocoronal Lya line, that 
swamps any intrinsic emission. 

Observations of similar dwarf galaxies show the 
deep Lya absorption seen in I Zw 18 is not unique. 
Two of the other most m etal-deficient galaxies k nown, 
SBS 0335-052 and Tol65 ( Thuan fc Izotov . 1997 ). have 
metallicities just a factor of 2 higher than iZwl8. 
Both also show broad Lya absorption profiles with 
equivalent widths of —20 to —30 A, no hints of Lya 
emission, and clear damping wings that imply Hi 
column densities above 2 x 10 21 cm -2 . Similar deep 
absorption is visible in the COS spect rum of low- 
metallicity dwarf SBS 1415+437 ( James et all 12014 . 
12 + log(0/H) a: 7.6), which also shows a similarly red- 
shifted bump of Lya in emission. These galaxies are very 
rare in the local universe, and one may question whether 
far-reaching conclusions may be drawn from them. Nev¬ 
ertheless, such objects may become more abundant at 
higher redshifts, and any complete theoretical picture 
of Lya must also include them. 


4.2 Galaxy Winds and Atomic Gas 
Kinematics 

In contrast to the strongly absorbing dwarfs the second 
BCG with a published HST Lya spectrum, Haro 2, was 
found to emit a strong Lya line with WL ya ~ 7 A (emis¬ 
sion part only, iLeoueux et all Il995l right panel of Fig¬ 
ure 0). This is particularly curious because Haro 2 is an 
order of magnitude more dust- and metal-rich than the 
dwarf galaxies. Comparing with the absorbing BCDs 
discussed in Section 14.11 the total column density of 
hydrogen along the line-of-sight is roughly the same, 
but with two important differences. Firstly, much less 
of that hydrogen column density is contributed by the 
neutral phase, although the measured column of Am 
> 10 20 cm -2 would still be sufficient to cause a damped 
absorption. Secondly, the absorption centroid of Lya is 
blue-shifted by ~ 200 kms^ 1 relative to the systemic 
frame of rest of the galaxy. The result of this first re¬ 
solved observation of a Lya emission line is an asymmet¬ 
ric profile that comprises a blue absorption component 
and a red emission peak, similar to the P Cygni profile. 
Even though the wing of the absorption profile is quite 
extended, the velocity offset is sufficient to shift the neu¬ 
tral medium partially out of resonance with Lya, and 
enable some of the Lya radiation to escape. 

The LIS lines intrinsic to Haro 2 were also found to 
be blue-shifted with respect to the systemic velocity 
and by the same velocity as measured from the Lya 
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absorption (« 200 kms : ). The peak of the Lya emis- regulation easily has strengthened. For gala xies in the 


( 201 lh P Cygni emission is ubiquitous while for the Ha- 
galaxies, redshifted Lya peaks and LIS lines blueshifted 
by around 100 kms' 1 are found among the emitters, 
while absorption lines (including Lya) consistent with 
zero velocity shift at 68% confidence are exhibited by 
the absorbers. 

Similar results are found among the Lyman 
alp ha Reference Sa mple (LARS; Section 15.21 


Rivera-Thorsen et al.1 2015]), which 


sion, however, is instead redshifted by 350 kms -1 , or ‘ Lyma n-break analog’ (LBA) samples of Heckman et al 
roughly twice the blueshift of the neutral medium. This 
led Lequeux et ah (1199511 to suggest that much of the 
Lya is able to avoid Hi absorption in Haro 2 because 
it does not see the atomic gas as static, and the red¬ 
shifted Lya emission supports a picture in which the 
Lya that is emit ted is ‘backscattered’ from a receding 
shell of Hi gas (jVerhamme et al.. 120081 ). Furthermore in 
Haro 2, diffuse soft X-ray emission covers and extends 
beyond the UV-bright, star-forming regions, which is 
produced by the mechani cal energy r ele ased b y the star 
formation episode ( Otf-Floranes et all 201 2(1. This X- 
ray emission is spatially consistent with an extension of 
the Lya emission in the 2D spectral image. In contrast 
I Zw 18, which shows only Lya absorption and static 
low-ionizatio n abso rption lines, is undetected at soft X- 
ray energies ( Ott et all 2005 ). The conventional picture 
for galaxy outflows is that cold gas is acc ele rate d by ex¬ 


panding hotter gas (e.g. Strickland et al. . 20041) . which 


supports scenario where feedback-driven outflows pro¬ 
mote the emission of Lya. 

This early picture easily g eneralized in a sa mple of 
eight local BCGs observed bv lKunth et al.l ( 19981 ). four 
of which show net Lya in emission and four absorp¬ 
tion. For those with net absorption, their OiA1302 and 
SinA1304 absorption lines lie close to the systemic ve¬ 
locity (within 25 kms -1 ), while the other four show 
outflowing gas with centroid velocities sh ifted by 60- 
180 kms -1 (see also iLeitherer et all 2013 ). This corre¬ 
lation does not necessarily imply a causal relationship 
and the fact that Lya is seen to be locally emitted where 
outflows are strong could also be explained also by ori¬ 
entation: transport models show more Lya to be emit¬ 
ted perpendicul ar to galaxy disks purely because of the 


Hi distribution (Vcrh amme et al.l . l2012t lLaursen et al 


2013 ), and winds are also stronger in the polar d irectio n 
because the pressure is also low er ie.g. iBland fc Tuliv . 
1988; Veilleux fc Ruukel . 2002 ). More compelling evi¬ 
dence for a causal association comes from the fact that 
all the Lya emission lines in the sample show P Cygni- 
like asymmetric profiles, indicating that photons are in¬ 
teracting directly in the outflowing medium. For this, it 
is much harder t o argue f or a non-causal relation. 

Galaxies in the iKunth et al. ( 1998 ft sample were orig¬ 
inally chosen to span a range of metallicities and dust 
contents, but both of these quantities were found to be 
secondary in governing Lya emission/absorption when 
compared to the presence/absence of outflowing neutral 
gas. 

COS Lya observations of local galaxies are ongoing, 
but already the instrument has far outdone the GHRS 
in terms of n umbers. Using l arger samples of both 
FUV-sel ected ( HeckmaneLaL, 20 111) and Ha-selected 
objects ( Wofford et all 20131 ) this picture of kinematic 


are summa¬ 
rized in the lower panels of Figure 0 The left-most 
panel shows an example where the ISM in front of the 
brightest nuclear star cluster (where the COS aperture 
is positioned) is static, and from where broad damped 
Lya absorption is also observed. From other regions of 
the galaxy, however, Lya emission is recovered and the 
galaxy becomes a weak Lya emitter (Wi, ya ~ 10A) in 
apertures that encompass the galaxy. The central panel 
shows an example where the neutral ISM is outflowing 
along the line-of-sight by around 250 kms -1 , the Hi 
absorption is similarly blueshifted, and a weak Lya 
emission feature is visible within the pointing of the 
COS. The right panel instead shows a galaxy where 
the atomic gas is outflowing at higher velocity still, 
and a very bright Lya emission line is visible with 
Wl ya ~ 40 A (80 A when including extended emission). 

Obviously if all neutral gas can be cleared from 
zero velocity, then Lya should escape unhindered. For 
the column densities probed using the Sin lines, this 
appears to be the case in the right panel of Fig¬ 
ure 0 where only more tenuous gas appears to re¬ 
main at Av « 0 to absorb Lya. Yet removal of cold gas 
from Av ss 0 n ot a re quirement for Lya emission, and 
Rivera-Thorsen et al. ( 20151 ) present several examples 


of galaxies that show strong Sin absorption at Av = 0, 
but also significant Lya emission (their Figure 8). How¬ 
ever while there is clearly gas that does not have a ve¬ 
locity shift, the centroids of the absorption profiles are 
offset - usually by < —50 kms -1 - which demonstrates 
that there is fast-moving gas that can Doppler shift Lya 
out of resonance with the static material. Additional 
support for this comes from the fact that the peak of 
the Lya profile is redshifted in every case. 

Figure 0 summarizes the situation, by showing the 
average velocity shift of low-ionization absorption lines 
compared with Wl ya measured in the same aperture 
(diameter of 1.9 and 2.5 arcsec with GHRS and COS, 
respectively). Outflow velocities span a range from —250 
to +50 kms -1 (net inflow), with a median value of 
—50 kms -1 . No galaxy with low-ionization gas moving 
in the velocity range —50 to +50 kms -1 shows signifi¬ 
cant Lya emission: all but two of these objects show ei¬ 
ther net absorption or emit Lya with a total EW below 
3 A. However galaxies with faster outflowing gas show 
a wide range of WL ya , which reaches up to 40A, with 
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Figure 5.: Example images and Ly a spectra of local galaxies. Upper panels show colour composite images, encoding 
Ha in red, FUV continuum in green, and Lya in blue. Dotted white lines indicate the size and position of the COS 
aperture. Lower panels show the corresponding spectra around Lya (dark blue) and SinA1260A (red). Objects 
are selected to illustrate the variety of Lya and LIS profiles. Note the different velocity scales, which extend to 
±7000 kms -1 in the left -most frame but just ±1500 kms -1 in the other two. The left object, LARS09, shows a 
completely absorbed Sin profile that is centred at zero velocity. This absorbing gas removes 300 kms -1 (FWHM) 
in the Sin line but completely damps the Lya line, resulting in a FWHM that is 20 times broader. The centre plot, 
LARS08, shows a similarly broad and saturated Sin profile, but one that is offset in velocity by ~ —250 kms -1 . 
The blue wing of the Lya absorption profile shows a similar shape and is offset in a similar way to the left panel, 
but in this case a redshifted Lya line is emitted. The right panel shows a an example of a strong Lya-emitter, in 
which the Sin absorption is blueshifted by « 350 kms -1 , is far from saturated, and a very bright Lya emission line 
is seen. 


an average near 10A. We can state with confidence that 
at least on small scales Lya emission is correlated with 
feedback of sufficient magnitude, that acts to accelerate 
the highest density neutral ISM along the line-of-sight. 
Note, however, that this does not mean it is the case 
that all galaxies with a strong outflow are Lya emitters 
- clearly there are galaxies with Av « —150 kms -1 and 
Wi, ya below 2 A. Furthermore while the correlation ex¬ 
ists, the causal mechanism by which feedback affects 
the transport has not necessarily been established. For 
example, whether feedback is simply shifting the Lya 
out of resonance by scattering from bulk-flowing gas, 


or the instigation of fluid instabilities that disrupt the 
ISM. 


4.3 Dissecting the Neutral Medium in Detail 


In addition to gas kinematics, the plethora of reso¬ 
nance absorption lines in the UV also provide a proxy 
for the fractional covering of the cold Hi me dium 


(e.g. Savage fe Sembach . 19961 : Pettini et al. . 2002 1. The 
principle is simple: strong resonance transitions are as¬ 
sumed to be saturated at normal metallicities and col¬ 
umn densities, and thus if the absorbing line does not 
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Figure 6.: The relationship between W\^ ya and outflow 
velocity of the atomic gas. All measurements are made 
within spectroscopic apertures with HST/GHRS and 
COS. Where the Lya absorption component dominates 
over emission Wl yoi is set to zero. Av is calculated 
from the average of velocity centroids of all observed 
LIS lines compared with the intrinsic velocity of the 
nebulae measured from optical emission lines. A cut of 
Av = — 50 kms -1 divides the sample in half by Av, 
segregating galaxies with fast outflows from those with 
weakly outflowing, static, or inflowing neutral media. 
Grey points with errorbars represent the average and 
standard deviation of galaxies of these two sub-samples. 


drop to the level of zero intensity then the observation 
hints that there may be multiple clouds inside the spec¬ 
troscopic aperture that do not fully cover the stellar 
sources of continuum radiation that lie beneath. Thus 
if there are direct sightlines between the observer and 
the ionized regions, Lya may escape unimpeded, and 
importantly, without frequency shift. Indeed if scatter¬ 
ing could be completely mitigated and dust confined to 
the cold gas phase, /^“ should be at least 1 — f c , where 
/ c is the Hi covering fraction. 

These methods have been used to place indirect limits 
on t he escape of ionizing radiation from sta r burst galax¬ 
ies ( Grimes et al. . 20091: iHeckman et all 201ll ). and 
have recently been verified by direct o bservations in the 


ioniz ing continuum with HST/COS (IBorthakur et al 


2014]). Similar tests, verified ag ainst Lya em i ssion, have 


been conducted at high-z (e.g. Ijones et~aI1 . 1201311 and 
low- 2 : COS studies focussed on Lya recently been pre¬ 
sented. For the majority of local UV-selected galax¬ 
ies the depth of the normally-saturated Sin lines indi¬ 
cates a covering fraction close to unity. However there 


is some deviation from this: a weak trend is seen for 
galaxies with fl^' a above 0.1 to be drawn from systems 
with SinA1260 A absorptio n lines that are not satu¬ 
rated (Rivera-Thoreen et ajJ, 120150 . Very well exposed 
continuum observations are needed to solve for covering 
fraction, but solutions include the possibility of f c < 1 
for the thickest gas neutral columns in Lya-emitting 
galaxies but not in the case of absorbing systems. 

Many of the resonant UV transitions, including the 
Sin discussed above, have an associated fluorescent 
transition at longer wavelength, denoted with a * 
(e.g. C 11 *), that provide additional diagnostics of the 
atomic medium (Prochaska et al.. 2011: Rubin et al 


201 lb Ijaskot fc Oevl . 201 4 Scarlata fc Panagia . 20151) . 
Since the absorption lines are resonant, they may be 
partially filled by scattered radiation, very similarly 
to Lya. However, the fluorescent transition associated 
with each line has a roughly similar Einstein A co¬ 
efficient to the resonant de-excitation, which implies 
that roughly half of absorbed photons should be emit¬ 
ted in the longer wavelength * line at every scatter¬ 
ing. In a symmetric, energy-conserving system without 
losses, absorption along the sightline must be balanced 
b y isotropic fluo r escent emission. 

Ijaskot fc Oevl ( 2014h recently presented spectra of 
two particularly interesting bright Lya emitters, with 
Whya between 70 and 150 A. Both of these objects are 
well-detected in the stellar continuum, but absorption 
lines of CliA1334 and SiiiA1260A lines are barely visi¬ 
ble. However the fluorescent counterpart of each tran¬ 
sition is clearly seen in emission, suggesting that the 
ISM is indeed partly covered or shows a low Hi col¬ 
umn density in these galaxies. A similar spectrum is 
that_of LARS 14^ ill ustra t ed in the lower-right panel of 
[5] ( Rivera-Thorsen et all 2015h . which shows a bright 
Lya line with a blue peak, incomplete Sin absorption, 
and a fluorescent emission line (seen at relative velocity 
of +1200 kms -1 ). 

Further information may be inferred from the profiles 
of Lya. While Lya and SiiiA1260 have cross sections 
of the same order of magnitude and become optically 
think at similar column density, the metal abundances 
imply that Lya may be absorbed by gas that is not vis- 
ible to metal absorp tion. Particularly in the galaxies of 
Jaskot fc Oevl ( 20141) and LARS 14, the Lya lines do not 


resemble the strongly asymmetric absorption+emission 
of P Cygni profiles that occur at high column densities 
of completely covered gas. Instead they show double- 
peaked profiles with narrow absorption at Av = 0. This 
indicates that there must be absorbing Hi that is not 
Doppler shifted, but that is also not of sufficient column 
density and/or metallicity to be seen in metal absorp¬ 
tion lines. This implies IVhi between 10 15 and 10 18 cm -2 
for normal ranges of metal abundance. 
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5 FROM SMALL TO LARGE SCALES: 

COMPLETING THE QUANTITATIVE 

PICTURE WITH HST IMAGING 

Are HST spectroscopic studies are performed in small 
apertures or narrow slits. Even for the most distant ob¬ 
jects discussed so far, the 2.5 arcsec entrance window of 
COS corresponds to a physical size of just 6 kpc, and ob¬ 
viously the apertures will sample much smaller scales in 
more nearby galaxies (e.g. just 130 pc in I Zw 18). Thus 
while providing a very rich picture about Lya and the 
ISM, spectroscopic measurements are restricted to the 
chosen sightlines: as shown in Figure [3] these will sam¬ 
ple only a fraction of the galaxies. These small apertures 
are necessary to get the high spectral resolution, but to 
capture a representative fraction of the Lya they are 
likely far too small, particularly when considering that 
Lya can scatter. To get this, spectroscopic observations 
must be complimented with large-aperture imaging. 

Thick columns of neutral gas have been observed in 
most of the starbursts mentioned so far, and spectro¬ 
scopic data strongly suggest that scattering removes 
Lya from the line-of-sight. An important question be¬ 
comes whether the bulk of the absorbed Lya is truly 
absorbed, or simply scattered to larger radii from where 
it is subsequently emitted. Such information is vital for 
the comparison with high -2 data, where ordinarily a 
flux and EW measurement may be available but little 
more, and standard practice is to adopt small aperture 
appropriate for point sources (~ 1 — 2 "). 

Lya imaging at z = 0 became efficient with HST 
when the Advanced Camera for Surveys (ACS) was in¬ 
stalled. Although it was technically possible with earlier 
cameras, the Solar Blind Channel of ACS has a total 
throughput of a few percent, and made such observa¬ 
tions feasible for the first time. Not only did this bring 
about the first resolved information on Lya (aside from 
very extended high -2 Lya-blobs), but also an immedi¬ 
ate leap to an angular resolution of w 0.1 arcsec. How¬ 
ever galaxies still need to be sufficiently redshifted in 
order for Lya not to be absorbed my Milky Way Hi 
(cz > 2500 kms -1 is sufficient) so obvious targets such 
as very well-studied nearby systems such as M82 or M33 
still cannot be observed. 

The first Lya images revealed a range of morpholo¬ 
gies. Near the central starbursts, Lya is seen in both 
emission and absorption, where it may vary between 
the two on sub-kpc scales with little obvious dependence 
upon local properties such as age or reddening. For ex¬ 
ample the first starburst with a p hotometrical l y cal i- 
brated Lya image, ESO338-IG04 ( Haves et all 20051) . 
shows a lane of Lya absorption that runs approximately 
E-W, and loops around one side of the galaxy only. 
This is seen at no other wavelength. In other regions 
jets of brighter Lya fuzz are visible, again seeming un¬ 
correlated with obvious signs of Lya production such as 


Ha emission. Obviously spectroscopic results will be a 
strong function of aperture placement. Furthermore, as 
soon as ACS was turned towards local starbursts, Lya 
halos were discovered to surround the star bursting re¬ 
gions ( Atek et al. . 2008 : Ostlin et al.l . l2009h . I now pro¬ 
ceed to discuss the results of small-scale resolved anal¬ 
yses and extended Lya halos, beginning with a case 
study of one system. 

5.1 Resolved Analyses: A Case Study of 
Haro 11 

Local lu minous BCG Haro 11 is a Lyman break ana¬ 
log (e.g. Grimes et al.. l2007l) . and emits Lya with a to¬ 
tal EW of 15 A. As shown in Figure [3 it comprises 
three main star-forming knots (labeled in the FUV im¬ 
age), all of which are bright in the UV and Ha, but 


only one of these condensations emits Lya ([Haves et al 


120071) . Following this breakdown of the galaxy, knots 
A and B (the west-most two) are by far the brightest 
in Ha and must produce the bulk of the Lya radia¬ 
tion, but both absorb at Lya. In contrast it is only 
the single easterly knot C, which is the faintest of the 
three in Ha but brightest in the UV, that locally emits 
its Lya. Haro 11 also emits a halo of Lya emission, 
centred around knot C, which can most easily be ex¬ 
plained if Lya is re-radiated after scattering in the sur¬ 
rounding neutral gas. Diffuse emission also surrounds 
the two easterly knots, but at the positions of the clus¬ 
ters absorption outweighs the emission, giving a nega¬ 
tive overall flux. While the Lya surface brightness of 
the halo is low, it is also very much larger than the 
UV continuum-bright regions, and in tot al contr ibutes 
~ 90 % of the total Lya flux ( Haves et all 12007 ). Re¬ 
sults inferred from small-aperture spectroscopic obser¬ 
vations will be a strong function of both the size and 
placement of the aperture. 

In Haro 11, Lya produced in knots A and B may 
still be emitted, and all the observation can say is that 
more radiation from the stellar continuum is absorbed 
locally than the sum of directly emitted Lya and any 
Lya scattered into the line-of-sight. Thus for a given 
pixel we still cannot say whether Lya is scattered and 
absorbed by dust locally, or whether it propagates some 
kpc and contributes to the halo emission. Remarkably 
in this three-region decomposition of Haro 11, is that 
the strongly emitting knot C is also the dustiest, and 
shows Eb~v ~ 0.4 magnit udes, w hile knot A in particu¬ 
lar is far less extinguished ( Atek et aLl . l2008l ). Under the 
simplistic assumption that dust plays a dominant role 
in regulating Lya visibility, this would be unexpected, 
although results from these Lya absorbing knots are 
reminiscent of the dwarf galaxies discussed in Section [4] 

Similar phenomena were noted throughout the sam¬ 
ple of nearby galaxies first observed with ACS. Specifi¬ 
cally ESO 338-IG04 was also found to exhibit a diffuse 
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Figure 7.: HST/ACS imaging of luminous blue compact galaxy Haro 11 ( Kunth et all 2003 : Haves et al. . 2007 1. 


Upper Left: far UV continuum; Upper Right: continuum-subtracted Ha; Lower Left: continuum-subtracted Lya. 
The Lower Right panel shows a composite of the other three images, with Ha in red, FUV in green, and Lya in blue. 
The physical scale in kpc is shown on the axes. In the FUV im age, the three main star-forming condensations (knots 
A, B, and C) are labeled in the nomencla ture o f Vader et al. (1199311 : the Lya-emitting clump with no detection at 
any other wavelength (Knot D in Kunth et al. , 2003h is labeled in the Lya image. 


Lya halo that dominates the Lya output (Haves et al 


200') ). and in a small sample of six local starbursts 
global Lya emission is invariably associa t ed with large 
scale halo emission ( Atek et al. . 2008t Ostlin et al 
l2009h . 


Nebulae produce Lya radiation intrinsically at an 
intensity of 8.7 times that of Ha (Section [2]), but when 
we contrast the local surface brightnesses of the two 
lines a very wide range of line ratios is found. These 
are illustrated in Figure 0 which contrasts the Ha and 
Lya surface brightness in Haro 11 pixel-by-pixel. We 
now discuss various regions of the diagram: 

• a. Lya can be emitted with fluxes similar to those 
expected from recombination. In Figure 0 above 
log Lya surface brightness of « —12.8, a region of 


proportionality is seen between Lya and Ha that falls 
right on top of the expectation value for Case B. These 
pixels correspond to knot C in Figure 0 Since this 
region of points is well defined and rather narrow, it 
can most easily be understood as Lya photons leaving 
the galaxy with little interaction with surrounding Hi. 
Lya photons emitted in these regions were most likely 
produced here. 

• b. Lya can be emitted with fluxes below those 
expected from recombination. The preponderance of 
points in Figure 0 lie below the case B line. This can 
be the result of two factors: Lya can be absorbed by 
dust, decreasing the Lya/Ha ratio (just as the Balmer 
decrement increases with dust in nebulae), or Lya can 
be scattered out of the line-of-sight by Hi. I.e. dust 
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Figure 8.: Comparison of Ly a and Ha surface photome¬ 
try for Haro 11. The abscissa shows the logarithmic Ha 
surface brightness, measured in individual pixels in the 
HST images, while the ordinate axis shows the corre¬ 
sponding surface brightness in Lya. I.e. contours show 
the density of points when comparing pixel values in 
the frames from Figure [TJ The red line shows the case 
B recombination line ratio of Lya=8.7xHa; intrinsi¬ 
cally astrophysical nebulae produce Lya and Ha radi¬ 
ation that follow this line. The plot is logarithmic, so 
only Lya emitting regions can be visualized. Histograms 
above and to the right show the overall the distribution 
of light emitted as a function of surface brightness. 


and Hi scattering act to move points down from the 
red line. Recall, also, that in this logarithmic plot Lya 
absorption cannot be visualized and more pixels are to 
be found at negative values of Lya. 

• c. Lya can be emitted with fluxes above those 
expected from recombination. Toward the upper left 
region of the diagram, Lya/Ha exceeds the value of 8.7 
expected for Case B. In this example, some pixels are 
10 times brighter in Lya than expected. This happens 
only at lower Ha surface brightness, and here Lya 
is spatially redistributed from elsewhere and emitted 
after scattering in the neutral ISM, resulting in the 
halo phenomenon discussed above and in the following 
section. Thus some of the Lya that shows Lya/Ha < 
Case B (point b) must be scattered and not simply 
attenuated by dust. 


5.2 Extended Halos 

Early imaging observations showed that in galaxies with 
net Lya emission, the dominant fraction comes from a 
component of large-scale extended emission that sur¬ 
rounds the star-forming regions. Indeed upwards of 50% 


of the Lya is typically emitted in halo regions that 
exte nded at least 10 kpc from the UV -bright clus¬ 
ters ( Haves et all 120051 : lAtek et al. . 2008ll . and in some 
galaxies this is the only Lya that emerges. The first 
large-scale Lya imaging survey of local st arbursts - t he 
Lyma n alpha Reference Sample (LARS, Ostlin ct al 


l2014h - has shown that extended halos are near ubiq - 
uitous in Lya emitting galaxies ( Haves et all 2013 1. 
Moreover LARS has enabled the first systematic sur¬ 
vey of the sizes of these halos and the comparison with 
other wavelengths. 

Continuum-subtracted Lya images show a wide range 
of morphologies. Typically they do qualitatively resem¬ 
ble those of the UV and Ha, but are more extended, 
and envelop the galaxies. On average halos have twice 
the linear size in Lya that the galaxy d oes in either 
the UV stellar continuum or Ha (Figured Haves et al 


2013 11. These estimates are made using the depth and 


redshift independent Petrosian radii, which are found to 
be below 15 kpc in Lya for LARS galaxies, with a me¬ 
dian value of 5 kpc. Furthermore, the extension of the 
Lya surface (£Lya> the ratio of Lya radius to Ha ra¬ 
dius) is not an independent quantity, and is correlated 
with a number of measured properties: notably ^Lya is 
anti-correlated with dust abundance, as demonstrated 
by the lower panel of Figure [9] and Lya extension is 
found to be larger at lower metallicit y and stellar mass . 
Radi ative transfer simulations (using IVerhamme et, all 


20121 ) show that this effect cannot be reproduced simply 


varying the dust content and what gives r ise to these 
extended halos is currently unclear. [Pardv et all ( 2014) 
have shown that higher £ Lyct is produced by galaxies 
with narrower 21 cm line-widths, but not the total mass 
in Hi (which neglecting mergers correlates with the line- 
width), possibly indicating that lower mass galaxies 
with less complex large-scale morphologies are the ones 
in which Lya scatters to the largest relative distances. 


6 LYMAN ALPHA SURVEYS AT LOW 
REDSHIFT 

A common criticism of some local Lya studies is that 
the sample selection is different from that of LAEs 
and LBGs at z > 2. The first generations of small 
HST studies w ere a ssembled largely from IUE atlas 
(IKinnev et all 119931 ). which in turn were selected from 
older objective prism or BCG surveys; in contrast most 
high-^ objects are selected by Lyman break techniques 
or Lya detection in narrowband filters. Several re¬ 
cent studies have begun to rectify t his: firstl y using 
HST t he Lyman Break Ana log samples ( Heckman et all 
2005 : Hoones et all 120071 based upon UV luminosity 
and co mpactness), t he Ly man a lpha Reference Sample 
(LARS; lHaves et alJl2013 . 20l4 based upon UV lumi¬ 
nosity and Ha EW), and studies undertaken with the 
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Figure 10.: The Lya luminosities probe d by 
galaxy survey s. Reference coding : C10= ICowie et al 


(2011); W 14 =IWold et al 
(2014); 1108=] Hand i et al 


(2011); Bll=|Blanc^^il. OOllb C12= ICantalupo et al 


( 20141) : H14= Haves et al 
(2008); Gll=G uaita et al 


( 20121) : H10= lHaves et al.1 ( 2010h . Narrowband surveys 
(Gil, C12, H10) have been artificially randomized by 
A 2 = 0.05 to aid visibility. 


UV continuum slope, /3 


Figure 9.: Extension of Lya ha los i n the LARS sam¬ 
ple; figure is taken from Haves et al. 1 201&) . The upper 
panel shows the Petrosian radius measured in Lya plot¬ 
ted against that measured in Ha. The two galaxies that 
globally absorb Lya are set to negative size. The dashed 
line shows the 1-to-l line, and Lya is on average twice 
the linear size of Ha. The lower panel shows the relative 
extension of Lya compared to Ha, £,hya- This example 
shows how the extension is anti-correlated with the UV 
continuum slope, /?, which is commonly used as a proxy 
for the dust content in galaxies. 


GALEX Satellite (De harveng et al.l . 2008; C owie et al 


2010l based upon selection by Lya emission). 

The GALEX satellite has been a vital contributer 
to low -2 Lya astrophysics. As well as FUV and NUV 
imaging channels, GALEX also had the capability to 
perform slitless spectroscopy across the same FUV 
and NUV bandpasses, providing low-resolution spec¬ 
troscopy of Lya for objects at 0.19 — 0.45 and 
2 « 0.65 — 1.25. In turn, this enables us to really 
survey the low -2 unive r se for Lya- e mittin g galaxies 
( Deharveng et all 20081 Cowie et all 2010h . and red - 
shifts around 1 ( Barger et al. . 2012t Wbld et all 120141) . 
in a manner very similar to those employed at high- 
2 . In the FUV channels GALEX LAE surveys cap¬ 
ture objects with NUV apparent magnitudes down to 


21.8 (AB); at 2 = 0.3 this corresponds to a SFR of 3.6 
Mg yr -1 , assuming the continuum is unobscured. 

Figure [TO] shows the Lya luminosities probed by 
various studies out to 2 ~ 4. As shown by the cyan 
squares and red upward-facing triangles, UV-selected 
LARS galaxies and GALEX LAEs occupy a very similar 
range of Lya luminosities. Contrasting these luminosi¬ 
ties with some z > 2 surveys, significant overlap is seen 
with the deepest g r ound-based observations : surveys of 
lHaves et al.l ( 20101) . C antalu po et al.l ( 2012 . both nar¬ 
rowband), and Rauch et, al.l ( 2008 . blind long-slit spec¬ 
troscopy, without aperture correction) overlap the local 
HST and GALEX samples at above the 50% level. 

The 2 ~ 1 LAEs discovered by GALEX, however, are 
an order of magnitude more luminous than the z < 0.4 
sample, and show almost no overlap with the more local 
objects. This is entirely a matter of distance, where the 
same limiting flux corresponds to a factor of 20 differ¬ 
ence luminosity between 2 = 0.3 and 1. This difference 
in luminosity may introduce some biases when com¬ 
paring the luminosity distributions of the two GALEX 
samples, as we show in the coming Section. These 2 ~ 1 
LAEs, however, do span a similar luminosity range to 
some of the l arger v olu me, shallower Lya-surveys at 
2 > 3 (e.g. [Blanc et al. . 2011 ). 

In the following Sections I discuss the most im¬ 
portant results from the (z) = 0.3 and ~ 1 studies 
with GALEX, regarding the numbers/evolution of Lya- 
emitting galaxies (Section ICT) . and statistical studies of 
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their properties fSection l6.2l) ; in part of course the prop¬ 
erties of galaxies are part of that evolutionary process. 

6.1 Evolution of Lyman alpha Galaxies into 
the Nearby Universe 

IUE demonstrated that Lya emission is rare in the 
nearby universe but GALEX could determine how rare. 
The first Lya luminosity functions (LFs) at z « 0.3 
showed that Lya emitting galaxies have become both 
fainter and less abundant than at high- z. Figure HT1 
shows the LFs measure d at z ~ 0.3 and 1 ( Cowie et all 
201flt Wold et al. . 20l4 respect ively 1 ), together with r e¬ 


A large fraction of this evolution takes place in the 
4.3 Gy r that elap ses between redshifts of 1 and 0.3 


cent measurements for z = 2.1 (Ciardullo et al., 2012). 


It should be noted that, while these LFs are the best 
that can be done with GALEX, the samples are not 
large: 119 at 2 ~ 0.3 and 141 at z ~ 1. However, because 
of GALEX’s wide field-of-view and the large continuous 
wavelength range provided by slitless spectroscopy, the 
cosmic volumes probed are in fact rather large, and cov¬ 
ering several 10 6 comoving Mpc 3 . Thus while the indi¬ 
vidual parameters in the Schechter function may not be 
very tightly constrained, the small number of galaxies 
is certainly because of the relative paucity of LAEs at 
lower redshifts. 

As discussed in Section the two GALEX samples 
cover different luminosity ranges (Figure |T0| . Moreover 
the dynamic range of the surveys is not large, and even 
for the less luminous 2 < 0.4 sample it is not possible 
to calculate the faint-end-slope (a) of the LF. Thus the 
LFs presented in Figure ITTI assume a, basing the as¬ 
sumption on measured values at z = 2 — 3 for the z ~ 1 
sample, and the a for local Ha-emitters for the z < 0.4 
sample. Note also that the faint-end slope of the Lya 
LF is also not very tightl y constr ai ned a t high redshift 
(reasons outlined in Dressier et all I2014I) . 

At z w 0.3 and to the UV limits mentioned in Sec¬ 
tion 1 the shape of the Lya LF closely resembles that 
of Ha and H/3, although is lower in normalization: Lya- 
emitting galaxies make up about « 15 % of the local 


Ha-selected gal axy population (IDeharveng et all 2008 


Hu et all 20' )9h and 1/20 the FUV counts at the same z 
(ICowie et all . 1 20 1(1) . This equat es to a volume-aver aged 
escape fraction of below 1 % ( Haves et all l201lh . al¬ 
though this / e L y“ will be slightly higher when con¬ 
sidering possible emission from galaxies not formally 
classed as LAEs (i.e. galaxies that emit weaker Lya 
with 0 < H'Lya < 20A). At 0 ~ 2 the average /^ a is 
~ 5 %, so the average Lya output of the whole cos¬ 
mic volume decreases 5-fold. Note that this decrease in 
the emitted fraction of Lya happens on top of the de¬ 
crease in the cosmic star formation rate density, which 
also drops by a facto r of 5-1 0 over the same change in 
redshift (Madau & Dickinson, 120141) , implying the Lya 
luminosity density of the local universe is far below that 
of z > 2. 


4.a V. yr tnat elap ses between redsmits 01 I and U.d 
( Wold et all 20141 ). At the higher redshift of 1 luminous 
Lya emitters are certainly in place, with luminosities 
equivalent to L* at z ~ 2 — 3, and the evolution between 
z ~ 0.3 and 1 can be well d escrib ed by a sim ple factor 
of 10 increase of L * in the Schechteil 1 19761) function. 
Over this redshift range the space density of LAEs does 
not appear to change but the galaxies simply scale up in 
luminosity. Between z = 1 and 2 (3 Gyr), </>* increases 
by an order of magnitude to produce the LFs observed 
in ground-based surveys. Arrows in Figure |TT] show how 
the evolution of the LF manifests as a drop in density, 
followed by a dimming. Between the peak in the cosmic 
SFRD at z ~ 2 — 3 and 2=1, the universe first acts 
to turn off a fraction of the Lya emitters, whereafter 
between 2 = 1 and the nearby universe the abundance 
is constant but the galaxies get fainter in line with the 
evolution in both the UV and Ha. 

Interestingly, this evolution is not strongly reflected 
in the equivalent width distribution of LAEs, which 
do es not evo lve as dramatically. The EW distribution 
of Guaita et al. ( 201ll) at 2 ~ 2.1 agrees well with the 
distribution at 2 ~ 1 (shown in the right panel of Fig¬ 
ure HU) despite the fact that the overall Lya-emitting 
fract ion has decrease d by a factor near 5 ( Cowie et al.l 
20 id contrasted with IShaplev et ah 20031) . Thus what¬ 


ever process is turning Lya-emitters off it does not af¬ 
fect the shape of the remaining EW distribution (i.e. 
galaxies with W\ jya > 20A). Of course the equivalent 
width distribution of the overall population changes sig¬ 
nificantly, as many galaxies drop below the canonical 
20A limit; the higher EW tail of the distribution re¬ 
mains largely constant. Note that it is not necessarily 
fair to conclude strong evolution in the EW distribution 
to 2 ~ 0.3 from Figure [Til as the GALEX FUV obser¬ 
vations are first continuum selected. This will lead to a 
fraction of continuum-faint objects with high EWs be¬ 
ing missed, which will extend the tail of the distribution. 
However because the bulk of the luminosity comes from 
low EW galaxies, the LFs are will be largely unaffected. 

In light of the extended Lya halos discussed in Sec- 
tion !5.21 we may ask whether much Lya also extends be¬ 
yond the spectroscopic extraction apertures of GALEX. 
These spectra are extracted usi ng an optimal model o f 
the point spread function CP SF . Morrissey et all 120071) . 
which has FWHM of ss 5". At 2 = 0.3 (2 = 1) this aper¬ 
ture corresponds to a spatial scale of 22 (40) kpc, and 
thus one dimensional spectra are summed over scales 
that exceed this. At 2 ~ 0.3 the aperture is 5 times 
the median Lya Petrosian radius in LARS, and while 
we do not known precisely how the Lya profiles be¬ 
have at larger radii, it is likely that in most cases 
GAL EX captures t he m ajority of the total flux. At 
2 = 1 [Barger et al. ( 20121 ) discovered a giant Lya blob, 
extended over an 18" diameter, but also measure such 
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Figure 11 .: Lef t: The evolution of the Lya luminosity fun ctions with redsh ift. Turquoise line shows the z ~ 0.3 LF 
of Cowie et al. ( 2010h and purple shows the 2 « 1 LF of Wold et al ( 2014 1. measured by GALEX FUV and NUV, 
respectively. The dark red line shows the LF of Ciardullo et al.1 (|2012h . measured at z = 2.1 using narrowband 
filters. Significant evolution is seen between the distant universe and the local one; grey arrows show the suspected 
m odes of ev ol ution in both luminosity and density (see text for details). The dotted line shows the Ha LF at z = 0.2 
of Tresse fe Maddox ( 1998ll . Right: the cumulative Lya equivalent w idth distribut i on at the sa me redshifts; data 
are taken from the same surveys, with the z « 2.1 distributions from lGuaita et ah (2010, 2011). 


objects to be extremely rare (one in the whole volume). 
Moreover, aperture sizes (~ 5" at z = 0.3) are equiva¬ 
lent to a 3" diameter aperture at 2 = 2, so if halo ex¬ 
tension does significantly affect the recovered Lya flux, 
it is likely by a similar factor as in high -2 observations. 


6.2 The Properties of Nearby Lyman alpha 
Galaxies 

Unlike most high -2 studies, both the GALEX and HST 
samples are sufficiently close that many of their proper¬ 
ties may be systematically measured. I assemble some 
of the key data obtained from these telescopes in Fig¬ 
ure [T2l and in this Section discuss what we have learned 
about galaxies that emit, and do not emit, Lya. 


6.2.1 Stellar Populations 

Over the Lya luminosity ranges probed by low -2 
Lya observations (~ 10 42 erg s _1 ), known LAEs have 
typical stellar masses on the order of ~ 10 9 M 0 
( Finkelstein et all 2011 ; Haves et, all 12014 ). This is sig¬ 
nificantly higher than estimates made for most LAE 
samples at high- 2 , while the Lya luminosities are 
equivalent or somewhat lower. This difference in the 
M/L(Lya) reflects the fact that the average nearby 
galaxy is more evolved, and has had substantially more 
time to build stellar mass than a similarly selected 
galaxy in the more distant universe. 


This evolution is also reflected in recovered stellar 
ages. However caution is needed here: when estimated 
from SED fitting age is a luminosity-weighted aver¬ 
age, and wavelength-dependent function of the included 
bands and the range of permit ted star-formation his¬ 
tories. Finkelstein et al.l ( 201 1') find a wide range of 
SED-fit ages with a median of ~ 300 Myr. However 
w hen su btrac ting a n underlying population of old stars, 
Haves et al. ( 20141) determine recent star-formation to 


have occurred in short bursts with luminosity-weighted 
ages up to a few tens of Myr. This suggests that Lya 
emitting galaxies, like starbursts in general, have temp 
temporarily elevated SFRs. The total Lya luminos¬ 
ity shows no dependency upon stellar age, but both 
W-'Lya and the relative throughput (/<^“) are higher 
in younger galaxies. 

The above would suggest that Lya escape fractions 
and equivalent widths are higher in galaxies with higher 
Ha equivalent wid th, w hich ro ughly measures the spe¬ 
cific SFR fsSFRh lCowie et all ( 201 lh indeed show that 
the fraction of Lya-emitting galaxies increases when 
progressively higher Ha EW thresholds are applied, 
finding al most 60% LAE s in sub-samples of 1 Fh 0 > 
250A. Haves et al.l (|2014h similarly find that all their 
galaxies with fl* a above 10% have Ha EWs above 
400A. Interestingly, however this effect is not visible 
when comparing the Lya/Ha ratio (or Wj jya ) with Wu a 
fFigurcfl2l) for the objects with W / L ya > 20 A. The most 
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Lyct emitters only 



UV magnitude EW(Ha) [ A ] 
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Emitters vs Non.Em. 
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log ( [N II] A6584/Hct ) 


Figure 12.: Scatter plots to the Left show the observed global Lya/Ha ratio m easured by G ALEX (pink squares) 


and H S T (dark blue circl e s). Fo r the GALEX samples, d ata are assembled from lCowie et all (1201ll) .IScailabi^t_al 


(2009), Finkelste in et al l (l201ll). and lAtek et al. 

(2009, 20141) and" Haves et all (|2oT3 . 20141) . Note that because of the differing quantities plotted on the abscissa 


(120141) : for the HST samples data are taken from Ostlin et a] 


the figures do not necessarily include the same number of points. Furthermore, since these points are chosen to 
have Lya in emission, the non-emitting population cannot be visualized. Upper Left: the absolute UV magnitude 
in the AB system. Upper Right: the Ha e quivalent width. Lower Left: the Ha/H/3 ratio, where the black line shows 
the effect of dust attenuation ( Calzetti et al . 2000l law), assuming Lya/Ha=8.7 and Ha/H/3=2.86. Lower Right: 
the [Nn]A6584A/Ha ratio (=N2 index). Histograms to the Right compare the frequency of GALEX-LAEs (pink 
filled histograms) and non-Lya-emitting galaxies (green outlined histogram), as a function of Wh<* {upper) and N2 
index {lower). 


convincing results emerge when we compare the average 
properties of LAE and non-LAE samples. 

6.2.2 Galaxy Morphology 

In UV-selected local galaxies, the objects with the high¬ 
est WLya and ff^ a are found to be among the more 
compact ones. LARS galaxies tend to have particularly 
compact UV morphologies and LAEs are found among 
those with Petrosian radii of ~ 1 kpc on average, sim i- 
lar to higher redshift results of [Malhotra et al.J (2012). 
Larger galaxies all show lower W\, ya . 

Galaxies hosting Lya-emitting starbursts represent 
a mixture of various types, although still a mixture 
that is distinct from the UV-selected galaxy popula¬ 
tion in general. The LAEs contain a higher fraction 


of compact galaxies and merging systems. More curi¬ 
ously, an enhanced frac tion of the GALEX LAE disks 
appear to be face on ( Cowie et al. . 201dl) . This phe- 
nome non is also expected from radiative transfer mod¬ 


eling ( Verhamme et al. , 2012 ; Laursen et al. , 20131) , and 
suggests that orientation effects may hide some disks 
from Lya selection. Unfortunately current samples are 
not sufficiently large to test Lya emission as a function 
of inclination angle. 

6.2.3 Interstellar Dust and Metals 

6.2.3.1 Dust 
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LAEs have bluer UV-optical colours tha n non- 


enhtting galaxies of the same magnitude (Cowie et al. 


2010l 201 ill but there is significant overlap in colours be¬ 
tween the emitting and non-emitting subsamples. UV- 
selected samples similarly exhibit higher Lya escape 
fractions (, f}Z a > 10 %) whe re UV colours are bluer 
(/3 < —1.8; lHaves et al.1 20141 ). This should be in part 
a reflection of the stellar age effects discussed in Sec¬ 
tion [H23 but may also be due to dust extinction that 
reduces the Lya throughput: Lya/Ha ratios and escape 
fractions (independent of age) are both higher for bluer 
galaxies, not only equivalent widths. 

The average Lya/Ha ratio found for_the GAL EX- 
selected LAEs is slightly above 2 (I Atek et all 12009a ; 
Scarlata et all 120091 : ICowie et all 201 Oh . although per¬ 
haps surprisingly, plots comparing Lya and Ha flux di¬ 
rectly show no general covariance over more than 1 dex 
in each quantity. Examining Lya/Ha however, the ra¬ 
tio decreases significantly with increasing Ha/H/3, and 
the conclusion that dust reduces the transmitted Lya 
in the LAE samples is shared over many studies. Indeed 
as shown in the lower left panel of Figure [T2k . the anti¬ 
correlation between Lya/Ha and Ha/H/3 is one of the 
few trends that is significant over the dynamic range of 
today’s surveys. 

The comparison of these line ratios with extinction 
laws reveals several curious features. Firstly, at the low¬ 
est Ha/H/3 ratios, most of the galaxies lie below the pre¬ 
dicted curves, and even in Lya-selected samples only 
~ 25% of Lya photons escape when Ha/H/3 is in the 
range 2.8-3.2. These galaxies emit less Lya than dust 
attenuation would predict. However the same is not true 
for dustier galaxies: at Ha/H/3 > 4 the mean Lya/Ha 
ratio is ~ 1, even though from this Balmer decrement we 
would would expect 97% of the Lya radiation to be ab¬ 
sorbed. Thus the locus of points in the Lya/Ha-Ha/H/3 
plane shows a trend that is much flatter than known ex¬ 
tinction laws, and the normalized Lya escape fraction 
(measured ///./“ divided by that which is expected for 
the derived dust content) increases with attenuation. 
Above Eb~v ~ 0.3 Lya emission beco mes on average 


several times stronger than expected (IScarlata et al 
20091 : Atek et al. . 20l4 Haves et~all 20141) . 


The most basic example of a simple screen of dust 
that redden s the nebular lines is i ncompatible with 
observ ation. Atek et all (2009a ) and Finkelstein et al.l 
(2009) invoked the Neufeld ( 199ll) geometry to explain 
this apparent enhancement of Lya, in which dust is 
embedded within the Hi clumps of a multiphase ISM, 
and Lya scattering prevents photons from encountering 
dust. Radiative transport simulations show that this 
effective ‘boost’ of Lya is very difficult to reproduce 
without rather contr i ved combinations of parameters 
( Laursen et al. . 20131 : Duval et al. . 2014 ). and the pre¬ 
dicted increase of Lya EW with measured attenuation 
is not observed. 


Scarlata et al. (20091) argue instead for a scenario that 
requires no such preservational scattering, but is still 
built upon a clumpy dust distribution (which may any¬ 
way follow the cold gas). This model does not require 
clumps to act as mirrors to Lya and nor d oes it predict 
Lya EW to rise with Eb-v- Similarly I Atek et al.1 ( 2014 ) 
point out that a galaxy is likely made up of many Hil 
regions with a large variety of optical depths and since 
we see only down to an optical depth of 1 at each wave¬ 
length, observed radiation in each line (and continuum) 
comes from regions of different sizes. At the dusty end 
of the galaxy distribution these latter two scenarios do 
n ot req uire an y sca ttering at all, and indeed recently 
Martin et al.l (2015) have shown that significant Lya 
emission can be detected from ultraluminous infrared 
galaxies (ULIRGs), from which it most likely escapes 
through holes in the ISM. 

6.2.3.2 Metal Abundance 


Dust can absorb Lya and metals cannot. At the 
UV luminosities probed by GALEX and LARS, the 
LAEs exhibit metallicities that overlap with the UV- 
continuum selected galaxies, but extend down to lower 
metal abundance. LAEs are on average deficient in 
nebular oxygen by about 0.4 dex (Figure fl2b). and 
are drawn mainly from a sub-populations with metal- 
licities of 12 + log(0/H) < 8.2. Indeed a remarkable 
result from the GALEX studies is that the N2 in¬ 
dex (=log([Nil]/Ha) segregates Lya-emitters from non¬ 
emitters more cleanly than Ha/H/3. 

Part of this apparent preference for low metallicity 
galaxies comes from an enhanced fraction of more com¬ 
pact irregular galaxies (Section 16. 2. 2D , which are lower 
metallicity in g enera l . Co mparing Lya EW with age 
and metallicity, ICowie et al.l 1 2011 ) suggest LAEs are 
drawn from a stage in the evolutionary sequence dur¬ 
ing which metals build up and Lya EWs decrease, in 
a similar w ay to that which would be expected for Ha. 
In addition Haves et al. ( 20141) find that not only does 
EW decrease with age but also the Lya escape frac¬ 
tion, which has no dependence upon stellar evolutionary 
stage. This shows that not only do older stellar popula¬ 
tions produce less Lya, but also that their ISM become 
more opaque to Lya; this would also support hypoth¬ 
esis of dust buildup, although current samples do not 
have the statistical power to say whether this is purely 
a dust effect. 


6.2.4 Neutral Gas 

While dust absorption is ultimately the process that 
can expunge Lya radiation from a galaxy, Hi scattering 
determines the path length of Lya between the nebulae 
in which it is produced and eventual emission. Here we 
discuss observations that probe the Hi phase directly by 
21 cm emission, telling us about the total amount of Hi 
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available for scattering, and its large scale kinematics. 
For a discussion of the Hi properties measured on small 
scales by absorption line studies, see Section 0] 

Masses in Hi, and the large scale Hi envelope of in¬ 
dividual galaxies can only be measured by 21 cm ob¬ 
servations, which currently limits us to the nearby uni¬ 
verse. However to avoid Lya absorption by the Milky 
Way and to separate the line from geocoronal emission, 
the lower limit of a galaxy’s recession velocity is cz < 
2500 km s _1 . By this distance of « 35 Mpc, Hi 21 cm ob¬ 
servations are already challenging ( Pardv et al .,2014). 

Regarding the Lya-absorbing dwarf galaxies, I Zw 18 
is known to lie within a huge Hi envelope that is many 
times its optical si ze, with a central co lumn density of 
Nm > 10 21 cm -2 I van Zee et all 1998h and total mass 
of 7 x 10 7 Mq. Thus if all the Hi were static with respect 
to the Hii regions, Lya would see upwards of 10 7 optical 
depths at line centre. SBS 0335-052 exhibits very simi¬ 
lar properties in bot h centr al L ya a bsorption and large 
static Hi envelope (Thuan et ah, 1999I) . In cases like this 
the fate of the Lya radiation is unclear: does the Lya 
scatter so many times that it is eventually absorbed by 
the small amount of available dust, or does the Lya 
eventually leak out in a large, low surface-brightness 
halo that cannot easily be detected? Without very large 
aperture imaging observations that cannot be obtained 
with existing facilities, conclusive answers are difficult 
to provide. 

Sixteen m ore st arbursts (Muv = —17.5 to —21, from 
Os tlin et al. l2009l and the LARS sample), complete with 


Lya imaging and spectroscopy from HST, have been 
observed at A = 21 cm with the Green Bank Telescope 
(GBT) and/or Karl G. Jansky Very Lar ge Array (VLA, 
Cannon et ah, 2004 : Pardv et al. . l2014h . Hi masses are 


between 10 9 to several 10 10 M 0 , which is or 10-100 
times as much as the dwarf galaxies , but for the seven 
galaxies where Hi detections are resolved the central 
column densities are comparable. In these samples the 
Lya-emitting objects with /^“ above 10% all have 
Hi masses below 4 x 1O 9 M 0 . While galaxies above this 
mass emit only a small fraction of their Lya photons, 
several of the less Hi-massive galaxies also show /^“ 
below a few percent, or in absorption. The same result 
is seen with dynamical masses derived from Hi line- 
widths. While it may be hypothesized that less Hi would 
permit more Lya photons to escape, it cannot be ruled 
out that this anti-correlation between f^. a and Mri is 
not simply reflecting the fact that higher fl J / c a is also 
seen at lower stellar mass (Section 16.2.11) and is not 
causally connected with Hi. 

Spatial resolution in these 21 cm observations is 
currently still low, even with the VLA, but in ev¬ 
ery case the Hi is resolved, and extends over more 
than 10 Petrosian radii in the UV. Most galaxies 
show signs of tidal interactions, both in their 21 cm 
morphologies and line shapes. Perhaps most interest¬ 


ingly ESO 338-IG04, IRAS 08 339+6517 (ICannon et al... 
120041) and Arp 238 (LARS 03, Pardv et all 20l4 ) show 
long Hi tidal streams and debris trails: in the cases 
of ESO 338-IG04 and IRAS 08339+6517 these streams 


stretch over 50-100 kpc to companion objects, while in 
Arp 238 a huge (M ~ 2 x 10 9 Mq) Hi body has been 
ejected that contains no apparent companion galaxy at 
the level of SDSS imaging. 


7 SYNTHESIZING THE 

OBSERVATIONAL DATA: WHAT 
MAKES A LYMAN ALPHA EMITTER? 

As we have seen in Sections 01 [5J and 16.21 current data 
show that Wj jya and f^ a are influenced by a large 
number of physical properties. In recent years much at¬ 
tention has been devoted to determining their order 
of precedence: i.e. whether covering fraction is more 
important than kinematic properties, which in turn is 
more important than dust reddening. Mostly we have 
been driven to find a way to predict the emergent Lya 
flux, EW or flf c a from a given set of conditions. We 
have now assembled significant samples of low- 2 : ob¬ 
servations, selected by Lya, Ha, and UV continuum, 
and have studied various subsets of them between the 
X-ray and radio, measuring all the physical properties 
discussed in Section [6] By combining this literature we 
now have substantial power to determine how Lya is 
regulated and test for the primary effects in homoge¬ 
neously selected galaxies that do and do not emit Lya. 
Simultaneously we need to explain: 


• Why only « 5% of local UV-selected galaxies 
(down to the SFRs of ~ 4M 0 yr -1 , that are within 
reach of GALEX) show W Bya above 20 A. 

• Why / e L s r and E B -v are anti-correlated in Lya- 
selected galaxies, but many galaxies significantly 
outlie this relationship; that many seemingly dust- 
free galaxies absorb their Lya, while many dusty 
galaxies exhibit f^/ c a that appears too high for 
their extinction. 

• Why higher WL ya and f^/ c a are found among 
galaxies with lower stellar mass and metallicity, 
and younger stellar age. 

• Why higher W Lya and / e L s £ a are found among more 
compact galaxies and face-on spirals. 

• Why Lya emission, at least on small scales, is fre¬ 
quently associated with galaxy outflows, and the 
spectroscopic line profile is almost always asym¬ 
metric. 

• Why among our highest EW Lya-emitting galax¬ 
ies we often infer covering fractions below unity 
(measured from, C 11 and Sin) or low Hi column 
densities. 
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In doing so we must remain cognizant of the fact that 
data have been assembled from different selection func¬ 
tions, measurements have been made in apertures that 
probe a variety of physical sizes, and that large scale 
halo emission may affect some, but not all, of our re¬ 
sults. 

We first address stellar properties, as these dictate 
the intrinsic Wj jya but do not affect transmission (i.e. 
they have no direct influence over fg/ c a , and only set 
a limit on the emergent Wi, ya ). Since Lya is repro¬ 
cessed ionizing continuum, Wl y Q . will exceed 20A only 
during for ages below ss 6 Myr ( Leitherer et all llQfM 


and Section EH), assuming an SSP. The SSP assump¬ 
tion is probably over-simplistic for whole galaxies but 
it does serve to illustrate that the Lya-briglit period 
cannot be sustained if the SFR is declining, and that 
episodes must be young. The Lya production must be 
most strongly correlated with the evolutionary stage of 
the stars. 

However the observed W\, ya is not strongly correlated 
with instantaneous SFR, and neither are the high-EW 
galaxies the most luminous in Lya. If anything the re¬ 
verse is true, at least within the current samples, and 
over the UV and Lya luminosity ranges probed at low-z 
high-EW galaxies are among the less luminous. Indeed 
similar can be said for galaxy samples selec ted b y other 


that show not only a lack of nebular Lya but also the 
absorption of continuum photons in the Lya resonance. 

If this dust were distributed purely as a screen sur¬ 
rounding the star-forming regions, the apparent over- 
luminous Lya observed in some LAEs cannot be ex¬ 
plained without severely modifying extinction laws. 
Lya would have to see at least the expected attenu¬ 
ation, plus an excess of absorption because of scatter¬ 
ing. To explain this we must invoke geometrical effects. 
Evenly mixing dust into the Hu regions would pro¬ 
duce an effective attenuation of the form oc 1 /Tdust 

at high Tdust, and asymptotically sets a lower limit 
to Lya/Ha ratio of ~ 2 (without scattering). How¬ 
ever such a dust geometry cannot make high a Ha/H/? 
ratio, which saturates for optically thick nebulae at 
ss 4. What can explain the simultaneous high Lya/Ha 
and Ha/H/3 ratio s is a clumpy distrib ution of dust, 
as advocated by Scarlata et al. ( 2009h . This m odel, 
originally implemente d by Natta fe Panagial (1198411 and 
Caplan &; DeharveneJ ( 1986 1. assumes dust to be dis¬ 


tributed in dense clumps and the effective attenuation 
law changes significantly with the average number of 
clumps along the line-of-sight. E.g. again without scat¬ 
tering, varying Td us t within 10 clumps may produce 
Lya/Ha> 2 with Ha/H/3«i 7, and can explain the line 
ratios seen in the dustiest local LAEs. The fact that 


emiss ion lines such as Ha or [QiilTA5007 (e.g|Atek_et_aL, Lya is seen at all from local ULIRGs ( Martin et al 


20111 ): the requirement for high EW line emission ne¬ 
cessitates not only ongoing star formation, but at a 
given SFR the UV continuum must be faint enough 
for the EW to be high. This at least partially explains 
why local LAEs are drawn from galaxi es wit h higher 
Ha EWs (as seen by, e.g. Cowie et all l201lh . and at 
lower stellar age and mass. C onsidering the well-kno wn 


2015h also suggests that Lya must find paths of low 


mass-metallicity relationship (Tremonti et ajj, 120041 ). it 


equally well explains why LAEs should be more preva¬ 
lent at lower metallicities. Note, however, that accord¬ 
ing to current spectroscopic data, this relationship does 
not extend indefinitely to the lowest metallicity, gas-rich 
dwarf galaxies; a point to which we will return below. 

Lya transfer however is not affected by any of the 
above considerations: stellar age, mass, and nebular 
oxygen abundance are not properties that influence ra¬ 
diation, although they may correlate with quantities 
that do. What matters from this point is the prop¬ 
erties of dust and Hi. f^ a is anticorrelated with the 
Eb-Vi as would be expected, but the spread is very 
large. Moreover the slope of the f^/ c a -EB-v relation¬ 
ship in LAE sample s is flatter t ha n w ould be expected 
for pure extinction ( Atek et all l2014h . indicating that 
the role of dust is diminished. At the dust-free end the 
under-luminous Lya can be most easily be explained by 
the presence of Hi, as it may scatter Lya photons many 
times, and increase the probability of dust absorption. 
This is supported by, for example, small aperture spec¬ 
troscopic observations of almost dust-free dwarf galaxies 


dust optical depth. 

Then we need to explain why the local dwarf galax¬ 
ies (e.g. I Zw 18, SBS 0335-052) centrally absorb at Lya 
despite their very low dust contents. Indeed HST spec¬ 
troscopy shows some absorption component in every ob¬ 
ject observed, demonstrating that Hi has effects that 
range between s mall d ips and very broad, damped ab¬ 
sorption features. [^norio^^gieeUaL ( 19991) presented 
an evolutionary sequence for the Lya spectral profile 
expected from a star cluster (assuming an SSP). This 
model assumes that during the earliest stages of a clus¬ 
ter’s evolution (< 2 Myr) local Lya absorption would 
be expected because the surrounding medium is com¬ 
pletely static. Over the subsequent ~ 4 Myr mechani¬ 
cal feedback from O star winds and the first supernovae 
accelerate the Hi outwards, producing P Cygni-shaped 
Lya profiles. The application of this model to dwarf 
galaxies suggests that their star-formation episodes may 
be too young and that feedba ck has not yet had tim e 


to accelerate the neutral ISM (IMas-Hesse et ah, 20031) . 


This scenario is probably an over-simplistic repre¬ 
sentation of whole galaxies, where stars form in dif¬ 
ferent regions over extended timescales, but neverthe¬ 
less the qualitative arguments may be helpful in un¬ 
derstanding the influence of galaxy winds. All the ob¬ 
served line profiles, including those of the local dwarf 
galaxi es, can in princip l e be u nified within such a sce¬ 
nario ( Mas-Hesse et all 120031) . In the complete litera- 
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ture of low -2 galaxies only one example shows symmet¬ 
ric Lya emission wi thout an obvious absorp tion compo¬ 
nent (Tol 1214-277. lThuan fc Izotov . 1997 ). and as Fig¬ 


ure [6] shows, an outflow in the neutral medium appears 


Lya emission (Kunth et al. 

19981: Wofford et al. 

Rivera-Thorsen et al.. 2015 

)■ 


2013 


8 PERSPECTIVES AND FUTURE 

DEVELOPMENTS: THE MOST 

PRESSING QUESTIONS 

With GALEX and HST Lya observations now running 
to over hundred galaxies for each telescope, local sam¬ 
ples are substantial. Nevertheless, several major ques¬ 
tions still remain to be answered. This closing Section 
is dedicated to outlining a handful of the most press¬ 
ing questions that, while perhaps challenging, can be 
addressed within current samples, or with possible ex¬ 
tensions with present-day facilities. 

8.1 What is the Atomic Gas Distribution in 
Star-Forming Galaxies? 

The spatial distribution of emitted and absorbed Lya 
must reflect some set of properties of the Hi gas. In 
very nearby galaxies, 21 cm observations have already 
revealed the structure of the cold atomic gas where, 
for example, the H i Nearby Galaxy S urvey (THINGS , 
Walter ct al. 2008 1 and VLA-ANGST (i'Ott et alllilli) 
surveys have found an atomic medium that is largely 
inhomogeneous and clumpy. It is clear that if Lya pho¬ 
tons were injected in these galaxies they would experi¬ 
ence more scattering in some regions than others. Ex¬ 
amining at least these 21 cm observations, we see that 
the Hi does not much resemble the shells and slabs in 
which many radiation transfer calculations are done. In 
order to empirically determine the effects that Hi and 
its distribution have on Lya emission requires resolved 
observations of individual targets in both Hi and Lya. 

Currently we may contrast average Lya surface 
brightness profiles (e.g. Haves et al.|, 2014) with those of 
the intensity at 21 cm (e.g. Bigiel fc Blit; . 120121) . These 
samples show both Lya and 21 cm intensity profiles that 
are best fit with Sersic profiles that are close to exponen¬ 
tial (n ~ 1), while UV/optical wavelengths in the LARS 
sample show significantly higher Sersic indices (n > 3). 
Does the Lya surface brightness trace the gas column 
density? While curious, this may be entirely coinciden¬ 
tal as the sample selection is very different in the two 
cases. Such observations established in the same galax¬ 
ies would be enormously instructive in interpreting the 
Lya halos of both galaxies at low and high redshifts. 

It is unfortunate that resolutions attainable in 21 cm 
and in the UV/optical are not well matched; sam¬ 
pling small physical scales at 21 cm requires very lo¬ 


cal galaxies, while the redshift requirements to observe 
Lya imply targets must lie b eyond several t e ns of Mpc. 
Progress has been made by Cannon et al.l (2004) and 
Pardv et al.1 ( 2014 ), but with synthesized beam sizes of 
« 15 arcsec in the best case (usually much coarser), 
such Hi observations can place only two resolution el¬ 
ements inside the linear size of the ACS/SBC camera. 
However VLA in configurations B and A can provide 
resolutions down to around 4 and 1.3 arcsec, respec¬ 
tively, albeit with a substantial increase in observing 
time. Observational programmes at higher resolution 
(VLA C and B configurations) are ongoing, but it is 
clear that large steps forward may be taken if Hi ob¬ 
servations are pushed to the highest spatial resolutions, 
and such observations are now essential. 


8.2 How is the Lya Spectral Profile Built? 

The total integrated Lya profile of a galaxy is built 
by emission from different regions, likely with differ¬ 
ing kinematics, orientation, and with different contri¬ 
butions to the total Lya. The spectral profiles of Lya 
measured in small apertures tend to be P Cygni, with 
radiation absorbed from the blue and re-radiated in 
the red. In these small apertures, Lya escape fractions 
are also low. However imaging tells us that these pho¬ 
tons are often not expunged from the system, but scat¬ 
tered back into the line of sight at different position, 
where large-aperture photometry measures significantly 
higher //(/“. It is an unfortunately common feature of 
absorption spectroscopy that we rarely know the line- 
of-sight distance between the emitting sources and the 
absorbing gas. 

Ideally we would like to know where the frequency re¬ 
distribution of Lya occurs, and whether the halo emis¬ 
sion shares a similar profile to that of the more central 
regions where the starburst is located. For example if 
frequency redistribution occurs close to ionizing clusters 
and Lya is singly scattered at large radii, the profile 
may indeed be similar over large distances. However if 
photons also get caught for many scatterings in halo gas 
then a new kinematic structure may be encoded in the 
Lya, or double-peaked profiles may be seen. 

To resolve this we may take high resolution spectra of 
the diffuse Lya-emitting regions in galaxy halos. These 
observations would be best supported by aperture- 
matched optical spectroscopy of Ha or H/3, to tightly 
constrain the rest velocity distribution of Lya, and ide¬ 
ally also the highest possible resolution 21 cm observa¬ 
tions (previous section). The best tool for this would 
be HST/STIS with narrow slits, which could provide 
resolutions of around 20 kms -1 , but again at the cost 
of lon g integrations in faint regions. Mas-Hesse et al.1 
( 2003i ) have shown that in IRAS 08339+6517, the blue 
wing of the Lya profile sets on at a similar wavelength 
over at least 10 kpc, suggesting the star forming regions 
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are surrounded by a Hi medium that may be rather ho¬ 
mogeneous in both space and velocity. However 21 cm 
observations find gas at much larger radii, where Lya 
has not yet been spectroscopically detected; whether the 
profile bends smoothly and traces the edge of a bubble, 
and how this effect may vary in different galaxies are 
all currently unknown. 

8.3 Ionization State of the Interstellar 
Medium 


Section 16.21 discusses the effect of a large number of 
galaxy properties on the emission of Lya. Many have 
been studied over time, but largely overlooked has 
been the effect of ionization state of the ISM. Specif¬ 
ically ionizing radiation from stars may not only pro¬ 
duce recombination nebulae but also heat the diffuse 
warm medium. Thus as well as producing Lya, the fur¬ 
ther propagation of LyC radiation may increase the 
ionization levels of the diffuse gas, lowering the opti- 
cal depth of th e Hll r egions to Lyman radiation (e.g. 
Pellegrini et ah, [2012 1. Given that the ISM of galax¬ 
ies may be very inhomogeneous, and that Lya may be 
absorbed close the nebulae in which it forms, the prop¬ 
agation of ionizing radiation could potentially have a 
large impact upon the first stages of Lya transfer. 

The ionization parameter (the number of hydrogen- 
ionizing photons per atom) governs the excitation 
of the gas, which is usually quantified observation- 
ally by the excitation parameter (an emission line ra¬ 
tio that contrasts high and low ionization species). 
This is done most effectivel y by t a king ratios of p 2 


and p 3 ions: Zastrow et al.l ( 2011 . 2013h performed 


‘ionization parameter mapping’ of several local star- 
bursts using the [Sin]A9069/[Sli] A6716 ratio, to iden¬ 
tify highly ionized cones that could signpost LyC 
emission. Similar integrated measurements of ‘Green 
Pea’ galaxies found high [Om]A5007/[Oii]A3727 ra¬ 
tios ( Jaskot fo Oevl . l2013ll that implies high ionization 
states, and indeed these galaxies have been found to be 
high-EW Lya emitters ( Jaskot fe Oevl 2014 : Henry et 
al., in preparation). 

Currently we do not know how these highly ionized 
regions affect Lya on small scales in the ISM. Re¬ 
solved Lya imaging has now been obtained for over 
50 galaxies with ACS, but as yet the ionizatio n struc¬ 
t ure has only been mapped in one of them. iBik et all 
( 2015h used VLT/MUSE observations of [Sill] and [Om] 
to map the ionization parameter and gas kinematics 
in ESO 3 38--IG04, and compare with the Lya imaging 
from lHaves et al. ( 20051 ). Two outflowing and highly- 
ionized conical regions are revealed, that approximately 
align with the brightest Lya regions. This would be 
consistent with a scenario in which Lya may propa¬ 
gate with less scattering through these highly ionized 
regions, but the observations of one system cannot es¬ 


tablish a causal relation, and a large sample of similar 
observations needs to be obtained. Pertinent observa¬ 
tional questions include whether Lya emission is sys¬ 
tematically enhanced in regions of high [Sin] / [Sli], or 
whether chimneys through the ISM may feed brighter 
regions of Lya emission in extended halos. 

8.4 Can we Predict Lya Observables from 
Other Information? 

Given the importance of Lya observations at both low 
and high redshift (see Section fTTIT) . a vital question be¬ 
comes whether we can predict the Lya escape fraction, 
equivalent width, or line profile from a given set of quan¬ 
tities. I.e. if we are given a dust content, a characteristic 
velocity for outflowing atomic gas, etc, do our predic¬ 
tions for Lya emission/absorption match reality? 

Given the absence of strong correlations involving 
Lya (for example in Figure fl2l) . one may be tempted 
to conclude that our understanding is not this sophisti¬ 
cated. However if, on the other hand, the Lya that we 
observe is governed by geometrical considerations such 
as viewing angle, then averaged over many galaxies the 
answer may be more encouraging. For example it has 
been shown that transfer models inside homogeneous 
shells can reproduce very wide ranges of line profiles 
I Schaerer et al.l ( 2011 ). similar to those that are observed 
globally at high-z; furthermore when coupled with semi- 
analytical models of galaxy formation, such models are 
able to reproduce the broad features of the Lya LF over 
a wide range of redshift ( Garel et all 12012). This may 
imply that our more general picture of Lya transport 
is correct. 

It is not clear whether current samples are sufficiently 
large, or span a high enough dynamic range in luminos- 
ity/SFR for such an exploration. However in the coming 
years, larger databases of global properties will become 
available; GALEX LAEs will remain at around 100 ob¬ 
jects; HST spectroscopic samples already exceed this 
while imaging observations remain somewhat smaller. 
This should permit statistical studies studies using lin¬ 
ear discriminant analyses, to determine how combina¬ 
tions of properties produce the observed Lya character¬ 
istics. Many such possibilities may be envisaged if the 
signal is strong enough and the samples are sufficiently 
large. 

More empirically motivated simulations can be per¬ 
formed using the wealth of data available for low-z 
galaxies as input. In such an approach, nature sets up 
the ISM instead of computers. For example, transport 
calculations have been run in synth etic g al axies output 
by hydro dynami c al simulations (Laursen et al.l . l2009t 
Verhamme et al. . 2012 : Yaiima et al. . 20141) . although 


as yet there has been no attempt to construct realistic 
input conditions of a galaxy based upon observation. 
Such an experiment is not easy, especially without de- 
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tailed knowledge of the Hi, but as that becomes avail¬ 
able the more the study becomes a possibility. With 
maps of the ionizing stellar population, nebular gas and 
ionization structure, Hi distribution including large and 
small scale kinematics, it will become possible to gener¬ 
ate sets of model galaxies that are based upon real sys¬ 
tems for which Lya observations have been obtained. 
Transfer simulations in such models will then recover 
the Lya morphology and spectral profile, that can be 
tested against observation. 

9 ACKNOWLEDGEMENTS 

I acknowledge the support of the Swedish Research Coun¬ 
cil, Vetenskapsradet and the Swedish National Space Board 
(SNSB). This research has made use of the NASA/IPAC 
Extragalactic Database (NED) which is operated by the Jet 
Propulsion Laboratory, California Institute of Technology, 
under contract with the National Aeronautics and Space 
Administration. I would like to extend warm thanks to my 
friends and collaborators for comments and feedback on the 
manuscript: John Cannon, Mark Dijkstra, Daniel Kunth, 
Peter Laursen, J. Miguel Mas-Hesse, Jens Melinder, Hector 
Otf Floranes, Ivana Orlitova, Daniel Schaerer and Anne 
Verhamme. Lucia Guaita is thanked for making equivalent 
width tables available for Figure 1111 Further I would like 
to thank Goran Ostlin, Claudia Scarlata, Sebastiano Can- 
talupo, Len Cowie, and many others, especially those from 
the 2013 NORDITA Lya workshops, for valuable and stim¬ 
ulating discussions. I thank the anonymous referee for sev¬ 
eral careful readings of the manuscript, and providing many 
thoughtful comments that have greatly improved the con¬ 
tent. 


REFERENCES 

Atek, H., Kunth, D., Hayes, M., Ostlin, G., & Mas- 
Hesse, J. M. 2008, A&A, 488, 491 

Atek, H., Kunth, D., Schaerer, D., et al. 2009a, A&A, 
506, LI 

—. 2014, A&A, 561, A89 

Atek, H., Schaerer, D., & Kunth, D. 2009b, A&A, 502, 
791 

Atek, H., Siana, B., Scarlata, C., et al. 2011, ApJ, 743, 
121 

Bacon, R., Accardo, M., Adjali, L., et al. 2010, 
in Society of Photo-Optical Instrumentation Engi¬ 
neers (SPIE) Conference Series, Vol. 7735, Society 
of Photo-Optical Instrumentation Engineers (SPIE) 
Conference Series 

Barger, A. J., Cowie, L. L., & Wold, I. G. B. 2012, ApJ, 
749, 106 

Bigiel, F.,_& Blitz, L. 2012, ApJ, 756, 183 

Bik, A., Ostlin, G., Hayes, M., et al. 2015, A&A, 576, 
L13 

Blanc, G. A., Adams, J. J., Gebliardt, K., et al. 2011, 
ApJ, 736, 31 


Bland, J., & Tully, B. 1988, Nature, 334, 43 
Borthakur, S., Heckman, T. M., Leitherer, C., & 
Overzier, R. A. 2014, Science, 346, 216 
Calzetti, D., Armus, L., Bohlin, R. C., et al. 2000, ApJ, 
533, 682 

Calzetti, D., & Kinney, A. L. 1992, ApJ, 399, L39 
Calzetti, D., Kinney, A. L., & Storchi-Bergmann, T. 
1994, ApJ, 429, 582 

Cannon, J. M., Skillman, E. D., Kunth, D., et al. 2004, 
ApJ, 608, 768 

Cantalupo, S., Lilly, S. J., & Haehnelt, M. G. 2012, 
MNRAS, 425, 1992 

Caplan, J., & Deharveng, L. 1986, A&A, 155, 297 
Chariot, S., & Fall, S. M. 1993, ApJ, 415, 580 
Ciardullo, R., Gronwall, C., Wolf, C., et al. 2012, ApJ, 
744, 110 

Cowie, L. L., Barger, A. J., & Hu, E. M. 2010, ApJ, 
711, 928 

—. 2011, ApJ, 738, 136 

Curtis-Lake, E., McLure, R. J., Pearce, H. J., et al. 
2012, MNRAS, 422, 1425 

Danforth, C. W., Tilton, E. M., Shull, J. M., et al. 2014, 
ArXiv e-prints, arXiv:1402.2655 
Davidsen, A. F., Hartig, G. F., & Fastie, W. G. 1977, 
Nature, 269, 203 

Deharveng, J. M., Joubert, M., & Kunth, D. 1986, in 
Star-forming Dwarf Galaxies and Related Objects, 
ed. D. Kunth, T. X. Thuan, J. Tran Thanh Van, 
J. Lequeux, & J. Audouze, 431-437 
Deharveng, J.-M., Small, T., Barlow, T. A., et al. 2008, 
ApJ, 680, 1072 

Dekel, A., Birnboim, Y., Engel, G., et al. 2009, Nature, 
457, 451 

Dijkstra, M. 2014, PASA, 31, 40 
Dijkstra, M., & Jeeson-Daniel, A. 2013, MNRAS, 435, 
3333 

Dressier, A., Henry, A., Martin, C. L., et al. 2014, ArXiv 
e-prints, arXiv:1412.0655 

Duval, F., Schaerer, D., Ostlin, G., & Laursen, P. 2014, 
A&A, 562, A52 

Finkelstein, S. L., Cohen, S. H., Malhotra, S., & Rhoads, 
J. E. 2009, ApJ, 700, 276 

Finkelstein, S. L., Cohen, S. H., Moustakas, J., et al. 
2011, ApJ, 733, 117 

Garel, T., Blaizot, J., Guiderdoni, B., et al. 2012, MN¬ 
RAS, 422, 310 

Giavalisco, M., Koratkar, A., & Calzetti, D. 1996, ApJ, 
466, 831 

Grimes, J. P., Heckman, T., Strickland, D., et al. 2007, 
ApJ, 668, 891 

Grimes, J. P., Heckman, T., Aloisi, A., et al. 2009, 
ApJS, 181, 272 

Guaita, L., Gawiser, E., Padilla, N., et al. 2010, ApJ, 
714, 255 

Guaita, L., Acquaviva, V., Padilla, N., et al. 2011, ApJ, 
733, 114 


PASA (2015) 

doi: 10.1017/pas. 2015.2014 


26 


M. Hayes 


Hartmann, L. W., Huchra, J. P., & Geller, M. J. 1984, 
ApJ, 287, 487 

Hartmann, L. W., Huchra, J. P., Geller, M. J., O’Brien, 
P., & Wilson, R. 1988, ApJ, 326, 101 
Hashimoto, T., Ouchi, M., Shimasaku, K., et al. 2013, 
ApJ, 765, 70 

Hayes, M., Ostlin, G., Atek, H., et al. 2007, MNRAS, 
382, 1465 

Hayes, M., Ostlin, G., Mas-Hesse, J. M., et al. 2005, 
A&A, 438, 71 

Hayes, M., Schaerer, D., Ostlin, G., et al. 2011, ApJ, 
730, 8 

Hayes, M., Ostlin, G., Schaerer, D., et al. 2010, Nature, 
464, 562 

—. 2013, ApJ, 765, L27 

Hayes, M., Ostlin, G., Duval, F., et al. 2014, ApJ, 782, 

6 

Heckman, T. M., Hoopes, C. G., Seibert, M., et al. 2005, 
ApJ, 619, L35 

Heckman, T. M., Borthakur, S., Overzier, R., et al. 
2011, ApJ, 730, 5 

Hill, G. J., Gebhardt, K., Komatsu, E., et al. 2008, in 
Astronomical Society of the Pacific Conference Series, 
Vol. 399, Panoramic Views of Galaxy Formation and 
Evolution, ed. T. Kodarna, T. Yarnada, & K. Aoki, 
115 

Hoopes, C. G., Heckman, T. M., Salim, S., et al. 2007, 
ApJS, 173, 441 

Hu, E. M., Cowie, L. L., Kakazu, Y., & Barger, A. J. 
2009, ApJ, 698, 2014 

Hu, E. M., Cowie, L. L., & McMahon, R. G. 1998, ApJ, 
502, L99 

Hummer, D. G., & Storey, P. J. 1987, MNRAS, 224, 
801 

James, B. L., Aloisi, A., Heckman, T., Sohn, S. T., & 
Wolfe, M. A. 2014, ApJ, 795, 109 
Jaskot, A. E., & Oey, M. S. 2013, ApJ, 766, 91 
—. 2014, ApJ, 791, L19 

Jensen, H., Hayes, M., Iliev, I. T., et al. 2014, MNRAS, 
444, 2114 

Jones, T. A., Ellis, R. S., Sclienker, M. A., & Stark, 
D. P. 2013, ApJ, 779, 52 

Kalberla, P. M. W., Burton, W. B., Hartmann, D., et al. 

2005, A&A, 440, 775 
Kennicutt, Jr., R. C. 1983, ApJ, 272, 54 
Kennicutt, Jr., R. C., & Kent, S. M. 1983, AJ, 88, 1094 
Keres, D., Katz, N., Weinberg, D. H., & Dave, R. 2005, 
MNRAS, 363, 2 

Kinney, A. L., Bohlin, R. C., Calzetti, D., Panagia, N., 
& Wyse, R. F. G. 1993, ApJS, 86, 5 
Kunth, D., Leitherer, C., Mas-Hesse, J. M., Ostlin, G., 
& Petrosian, A. 2003, ApJ, 597, 263 
Kunth, D., Lequeux, J., Sargent, W. L. W., & Vialle- 
fond, F. 1994, A&A, 282, 709 
Kunth, D., Mas-Hesse, J. M., Terlevich, E., et al. 1998, 
A&A, 334, 11 


Lacy, J. H., Malkan, M., Becklin, E. E., et al. 1982, 
ApJ, 256, 75 

Lanzetta, K. M., Bowen, D. V., Tytler, D., & Webb, 
J. K. 1995, ApJ, 442, 538 

Laursen, P., Duval, F., & Ostlin, G. 2013, ApJ, 766, 124 
Laursen, P., Sommer-Larsen, J., & Andersen, A. C. 
2009, ApJ, 704, 1640 

Lebouteiller, V., Heap, S., Hubeny, I., & Kunth, D. 
2013, A&A, 553, A16 

Leitherer, C., Clrandar, R., Tremonti, C. A., Wofford, 
A., & Schaerer, D. 2013, ApJ, 772, 120 
Leitherer, C., Schaerer, D., Goldader, J. D., et al. 1999, 
ApJS, 123, 3 

Lequeux, J., Kunth, D., Mas-Hesse, J. M., & Sargent, 
W. L. W. 1995, A&A, 301, 18 
Madau, P., & Dickinson, M. 2014, ARA&A, 52, 415 
Malhotra, S., Rhoads, J. E., Finkelstein, S. L., et al. 
2012, ApJ, 750, L36 

Martin, C., Moore, A., Morrissey, P., et al. 2010, 
in Society of Photo-Optical Instrumentation Engi¬ 
neers (SPIE) Conference Series, Vol. 7735, Society 
of Photo-Optical Instrumentation Engineers (SPIE) 
Conference Series 

Martin, C. L., Dijkstra, M., Henry, A., et al. 2015, ApJ, 
803, 6 

Mas-Hesse, J. M., Kunth, D., Tenorio-Tagle, G., et al. 
2003, ApJ, 598, 858 

Meier, D. L., & Terlevich, R. 1981, ApJ, 246, L109 
Morrissey, P., Conrow, T., Barlow, T. A., et al. 2007, 
ApJS, 173, 682 

Natta, A., & Panagia, N. 1984, ApJ, 287, 228 
Neufeld, D. A. 1991, ApJ, 370, L85 
Oke, J. B., & Zimmerman, B. 1979, ApJ, 231, L13 
Osterbrock, D. E. 1962, ApJ, 135, 195 
Ostlin, G., Hayes, M., Kunth, D., et al. 2009, AJ, 138, 
923 

Ostlin, G., Hayes, M., Duval, F., et al. 2014, ApJ, 797, 

11 

Oti-Floranes, H., Mas-Hesse, J. M., Jimenez-Bailon, E., 
et al. 2012, A&A, 546, A65 

Ott, J., Walter, F., & Brinks, E. 2005, MNRAS, 358, 
1453 

Ott, J., Stilp, A. M., Warren, S. R., et al. 2012, AJ, 144, 
123 

Pardy, S. A., Cannon, J. M., Ostlin, G., et al. 2014, 
ApJ, 794, 101 

Partridge, R. B., & Peebles, P. J. E. 1967, ApJ, 147, 

868 

Pena-Guerrero, M. A., & Leitherer, C. 2013, AJ, 146, 
158 

Pellegrini, E. W., Oey, M. S., Winkler, P. F., et al. 2012, 
ApJ, 755, 40 

Pentericci, L., Vanzella, E., Fontana, A., et al. 2014, 
ApJ, 793, 113 

Pettini, M., Rix, S. A., Steidel, C. C., et al. 2002, ApJ, 
569, 742 


PASA (2015) 

doi: 10.1017/pas. 2015.2014 



Nearby Lyman alpha galaxies 


27 


Pritchet, C. J. 1994, PASP, 106, 1052 
Prochaska, J. X., Kasen, D., & Rubin, K. 2011, ApJ, 
734, 24 

Raiter, A., Schaerer, D., & Fosbury, R. A. E. 2010, 
A&A, 523, A64 

Rauch, M., Haelrnelt, M., Bunker, A., et al. 2008, ApJ, 
681, 856 

Rivera-Thorsen, T. E., Hayes, M., Ostlin, G., et al. 
2015, ApJ, 805, 14 

Rubin, K. H. R., Prochaska, J. X., Menard, B., et al. 
2011, ApJ, 728, 55 

Savage, B. D., & Sembach, K. R. 1996, ARA&A, 34, 
279 

Scarlata, C., & Panagia, N. 2015, ApJ, 801, 43 
Scarlata, C., Colbert, J., Teplitz, H. I., et al. 2009, ApJ, 
704, L98 

Schaerer, D. 2003, A&A, 397, 527 

Schaerer, D., Hayes, M., Verhannne, A., & Teyssier, R. 

2011, A&A, 531, A12 
Schechter, P. 1976, ApJ, 203, 297 
Shapley, A. E., Steidel, C. C., Pettini, M., & Adel- 
berger, K. L. 2003, ApJ, 588, 65 
Stark, D. P., Ellis, R. S., Chiu, K., Ouchi, M., & Bunker, 
A. 2010, MNRAS, 408, 1628 
Steidel, C. C., Bogosavljevic, M., Shapley, A. E., et al. 
2011, ApJ, 736, 160 

Strickland, D. K., Heckman, T. M., Colbert, E. J. M., 
Hoopes, C. G., & Weaver, K. A. 2004, ApJ, 606, 829 
Takada, M. 2010, in American Institute of Physics 
Conference Series, Vol. 1279, American Institute of 
Physics Conference Series, ed. N. Kawai & S. Na- 
gataki, 120-127 

Tenorio-Tagle, G., Silich, S. A., Kunth, D., Terlevich, 
E., & Terlevich, R. 1999, MNRAS, 309, 332 
Terlevich, E., Diaz, A. I., Terlevich, R., & Vargas, 
M. L. G. 1993, MNRAS, 260, 3 
Thuan, T. X., & Izotov, Y. I. 1997, ApJ, 489, 623 
Thuan, T. X., Lipovetsky, V. A., Martin, J.-M., & 
Pustilnik, S. A. 1999, A&AS, 139, 1 
Tremonti, C. A., Heckman, T. M., Kauffmann, G., et al. 
2004, ApJ, 613, 898 

Tresse, L., & Maddox, S. J. 1998, ApJ, 495, 691 
Tumlinson, J., Thom, C., Werk, J. K., et al. 2013, ApJ, 
777, 59 

Vader, J. P., Frogel, J. A., Terndrup, D. M., & Heisler, 
C. A. 1993, AJ, 106, 1743 
Valls-Gabaud, D. 1993, ApJ, 419, 7 
van Zee, L., Westpfahl, D., Haynes, M. P., & Salzer, 
J. J. 1998, AJ, 115, 1000 
Veilleux, S., & Rupke, D. S. 2002, ApJ, 565, L63 
Verhannne, A., Dubois, Y., Blaizot, J., et al. 2012, 
A&A, 546, Alll 

Verhannne, A., Schaerer, D., Atek, H., & Tapken, C. 
2008, A&A, 491, 89 

Verhannne, A., Schaerer, D., & Maselli, A. 2006, A&A, 
460, 397 


Walter, F., Brinks, E., de Blok, W. J. G., et al. 2008, 
AJ, 136, 2563 

Wofford, A., Leitherer, C., & Salzer, J. 2013, ApJ, 765, 

118 

Wold, I. G. B., Barger, A. J., & Cowie, L. L. 2014, ApJ, 
783, 119 

Wu, C.-C., Boggess, A., & Gull, T. R. 1980, ApJ, 242, 

14 

Yajima, H., Li, Y., Zhu, Q., et al. 2014, MNRAS, 440, 
776 

Zastrow, J., Oey, M. S., Veilleux, S., & McDonald, M. 
2013, ApJ, 779, 76 

Zastrow, J., Oey, M. S., Veilleux, S., McDonald, M., & 
Martin, C. L. 2011, ApJ, 741, L17 


PASA (2015) 

doi: 10.1017/pas. 2015.2014 



